Screening in Hot 
Non-Abelian Plasma 

(PhD Thesis) 

Peter Petreczky 
June 25, 1999 



Eotvos Lorand University 

Budapest 

Supervisor: 

Andras Patkos 



Contents 



1 Introduction 



2 Review of non-Abelian Screening 

2.1 Chromoelectric Screening at Leading Ordei 



2.1.1 Chromoelectric Screening in Linear Response Theory 

2.1.2 Debye screening from classical kinetic theory 



2.1.3 Polyakov Loop Correlator and the Debye Screening! 



2.2 Infrared divergencies Associated with Static non-Abelian Magnetic Field; 



2~3 Dimensional reduction 



2.4 The Non-Abelian Debye Screening beyond Leading Order 



2.4.1 Non-Abelian Debye Screening in One-Loop Rcsummed Perturbation 



Theory and the Magnetic Mass 



2.4.2 Non-Perturbative definitions of the Debye mass 



2.5 Screened Perturbation Theory 



3 Coupled Gap Equation for the Screening Masses in Hot SU(N) Gauge 
Theory 



6 

6 
6 
7 
8 
11 
12 
13 

13 
15 
16 



19 



3.1 Gauge non-invariant resummation schemej 19 

3.2 Gauge Invariant Approach! 22 



3.3 Contribution of non-Local Operators to the Gap 

Equation! 23 



3.4 What Have We Learnt from the Coupled Gap Equations ? 25 



4 Screening Masses of Hot SU(2) Gauge Theory from Monte-Carlo Simu- 



lation of the Lattice 3d Adjoint Higgs Model 



27 

1.1 The 3d SU(2) Adjoint Higgs Model on Lattice] 28 

4.2 Determination of the Propagators and the Screening Masses on Lattice] ... 29 

4.3 Numerical Results for the Propagators of the 3d SU(2) Higgs Model .... 31 

4.4 Magnetic Mass in 3d SU(2) Pure Gauge Theory] 34 



4.5 Gauge Invariant Correlators and the Additive Constituent Gluon Mod"e]| ... 35 

4.6 Gauge Dependence of the Screening Masses! 37 



5 Coupled Gap Equation for SU(2) Higgs Model 




5.1 The Extended Gap Equations 






5.2 Numerical Results 




5.3 Screening Masses in the Symmetric Phase with a Gauge Invariant Resumma- 


tion Scheme 



6 Conclusions 



A Calculation of Different Diagrams in SU(JN) Adjoint Higgs Model 



40 

41 
44 

48 
51 
53 



[B 3d Adjoint Higgs Model on Lattice 



57 

B.l Phase Diagram of 3d Adjoint Higgs Model and Dimensional Reduction! ... 57 
B.2 Extracting the Screening Masses from the Propagators] 61 



B.3 The Dependence of the Screening Masses on the Choice of the Temperature 

Scale and Lattice Spacing! 62 



C Self-energy Contributions in 3d fundamental + adjoint Higgs model 



64 



3 



1 Introduction 



According to the standard cosmological model after the Big-Bang the density and the temper- 
ature of the early Universe were very high and it consisted actually of a plasma of different 
elementary particles. This fact is a great opportunity for particle physics since the early 
Universe thus serves as a laboratory for particle physics. In order to describe the matter 
under such extreme conditions a new formalism has been developed in the past 25 years 
which is nowadays called finite temperature quantum field theory (see Refs. |TJ, Q for detailed 
monographies on the subject). 

As the expanding Universe was cooling down it underwent several phase transitions. Two 
of these phase transitions, namely the electroweak phase transition (EWPT), which occured 
when the temperature was T ~ 200GeV and the QCD deconfinement phase transition, which 
happened at temperature T ~ 200MeV r are of special interest. The underlying theories are 
well established theoretically and tested experimentally []. The neighborhood of EWPT 
is relevant for the understanding the evolution of the baryon asymmetry of the Universe 
(see 0, U for reviews). The theoretical study of QCD deconfinement phase transition is of 
great practical interest because there is a hope to create quark-gluon plasma in relativistic 
heavy ion collisions. Actually, the possibility of producing the quark-gluon plasma state 
experimentally was one of the main motivations for detailed studies of this phase transition. 
One may hope to suggest signals of its existence which can be tested experimentally and to 
measure some of the thermodynamical properties of the new phase which can be confronted 
with the expectations based on theoretical models. 

One of the most natural signatures of the deconfined phase is the screening of static 
chromoelectric fields. Actually this was the first signature based on which the existence 
of the new phase was proven theoretically H |(|. More precisely in the confined phase the 
potential between heavy quark-antiquark pair rises linearly with the distance while in the 
deconfined phase this potential is exponentially screened, with an inverse screening length 
equal to the Debye screening mass. The screening of static chromoelectric fields is also of 
great phenomenological interest in detecting the quark-gluon plasma J7^ ||. 

In the ideal case when the temperature is sufficiently high the perturbative approach 
seems to be adequate. The coupling constant is small as the consequence of asymptotic 
freedom : at high temperature the only scale which is available is the temperature scale 
itself. The coupling constant is small at large temperature T: g{T) ~ 1/ ln(T/ Aqcd)- Thus 
naively one would think that the quark-gluon plasma is close to the ideal gas of quarks 
and gluons (the situation is similar for the high temperature electroweak plasma). However, 
this simple picture is spoiled by infrared divergencies. It turns out that beyond the natural 
length scale 1/T there are different screening legth scales, namely the electric, or the Debye 
scale 1/gT and the magnetic length scale l/g 2 T. While the electric scale is similar to some 
extent to the Debye screening scale of QED plasma the magnetic scale is present only in 
the non-Abelian case. Over the past few years important progress has been made in our 
understanding of the role of these scales 0, |H], [13], Electric screening scale is 

very important in phenomenological investigation of the thermodynamics of the quark-gluon 
plasma and calculating the signals of the existance of this state [15|, |16[ . The existence of 

1 The Higgs particle is not yet discovered and supersymmetric extensions of the electroweak theory might 
be relevant 
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the electric screening scale sets the boundary of applicability of the of naive perturbative 
approach and provides a possibility to cure some infrared divergencies and to perform a 
consistent resummation of the perturbative series [|| |T3. In the non-Abelian gauge theory 
there are, however, also infrared divergencies due to the presence of the magnetic mass scale. 
This leads to the breakdown of the perurbative approach beyond some order of perturbation 
theory [T]J . It is believed that the existence of the magnetic mass scale is responsible for the 
termination of the I s * order phase transition in the electroweak theory. This occurs when 
the Higgs mass is larger than some critical value |18| . 

In the calculation of static properties at high temperature the most efficient way to deal 
with the presence of different legth scales is the procedure of gradual reduction of the theory 



22|, |23|, |4[ |. In the past few years considerable progress has been made in understanding 



the electroweak phase transition based on this approach pi O, 12, 19, BTJ (see Ref.fi] for a 



review). Dimensional reduction was also successfully applied to the deconfined phase of QCD 



25| , f2q , |27| , |28|| . The dimensionally reduced QCD was used for non-perturbative definition 
of the chromoelectric screening p9[ |30| . However, more recent investigations revealed some 
problems in application of the dimensional reduction in QCD [Q. 

The present thesis is the comprehensive summary of my work devoted to the investigation 
of field theoretical screening phenomena by resummation of perturbative series and lattice 
Monte-Carlo technique. Also I have examined the consistency and precision of dimensional 
reduction. The organisation of the presented thesis is the following. In section 2 I am 
going to review the present status of chromoelectric and chromomagnetic screening and 
their applications to high temperature non-Abelian theories. Section 3 devoted to the self- 
consistent determination of the screening masses of hot SU(N) gauge theories. In section 
4 the validity of dimensional reduction is carefully examined and the screening masses are 
defined by non-perturbative (lattice Monte-Carlo) technique. In section 5 the screening 
masses of SU(2) Higgs model will be examined and a new procedure for improving upon the 
standard of dimensional reduction will be proposed. 

Sections 3, 4 and 5 contain the results obtained by the author in collaboration with F. 
Karsch, M. Oevers, A. Patkos and Zs. Szep. 



2 The gradual reduction approach is based on integrating out consecutively the heavy fields in the Eu- 
clidean path integral. In the first step the non-zero Matsubara modes with typical mass ~ nT are integrated 
out leaving an effective 3d theory. Therefore this approach is also referred to as dimensional reduction. 
For small enough couplings futher reduction is possible by integrating out the static Aq field (temporal 
component of the gauge field) with a mass ~ gT. 
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2 Review of non-Abelian Screening 



In this section a detailed review of non-Abelian screening phenomenon will be given. The 
first subsection deals with definition of non-Abelian Debye screening at leading order of per- 
turbation theory. Also a nice physical picture for the non-Abelian chromoelectric screening 
will be given there. In the following subsections we will discuss infrared problems of finite 
temperature gauge theories and will introduce the idea of dimensional reduction. Finally, 
the problem of chromoelectric screening beyond the leading order of perturbation theory will 
be reviewed. 



2.1 Chromoelectric Screening at Leading Order 
2.1.1 Chromoelectric Screening in Linear Response Theory 

The phenomenon of chromoelectric screening could be easily understood as the field induced 
by static colour source placed into the plasma. First let us study the response of the plasma 
to an arbitrary weak external field. The classical source which induces this field is j d . The 
perturbation induced by this source is 

V = J d 3 xf cl A,{x) (1) 

The classical field induced by this extremal source is 

A,(x) =< A^x) >= -i J d 3 x'D«(x - x')f d (x), (2) 

where D^ u (x — x') is the retarded gluon propagator. In the momentum representation the 
retarded gluon propagator could be written as |2[ 

B R = - P T + - P L jt K » Kv (%\ 

K 2 - U T {K) ^ K 2 — U L (K) ^ K 4 ' y ' 



where the following projectors were introduced 



P^(K) = - k -§-)5l (4) 



K 2 

In the case of static chromoelectric source J^(x) = Q5q5 3 (x) the induced potential is 



P^(K) = (<V - -^f) - Pi. (5) 



/d 3 k e ikx 
M ^ + n £ fo = o,k) = (6) 



(27r) 3 A; 2 + n L (A; = 0,k) 
Q f°° e ikr - e- ikr kdk 



2-n 2 2ir k 2 + U L (0,k) 



(7) 
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At leading order in the high temperature limit 

U L (k = 0, k) = U 00 (k = 0, k) = i(iV + N f /2)g 2 T 2 = m 2 m . (8) 

for SU(N) and Nf fermion flavour. Now Eq. (||) could be easily eveluated by closing the 
contour in the upper and lower half complex k-plane and picking up contribution from simple 
poles at k — iimco, which yields 

Airr 

Thus we can see that the field induced by static colour charge in the plasma is described by 
screend Coulomb potential. This phenomenon is called chromo electric or Debye screening ^ 
and the corresponding inverse screening length is called the Debye mass. 



2.1.2 Debye screening from classical kinetic theory 

The leading order Debye screening mass is determined by the high temperature limit of the 
one-loop self energy diagram. It was shown that the contribution of such diagrams, the 
so-called hard thermal loop can be described by an effective kinetic theory f35 |. Moreover 



a simple classical kinetic description could be given for the associated phenomena [36]. The 
idea of the classical kinetic description is the following. 

The large momentum modes (p > T) of the gauge and fermion fields are described as 
classical particles with some one-particle distributions. If there are no external fields all 
particles have equilibrium distributions. If one applies a weak external field to the plasma, 
the one-particle distributions change and induced currents appear which are proportional to 
the change of the on-particle distributions. The induced currents create some classical field 
which is the response of the plasma to external perturbation. This idea should be worked 
out self-consistently. For this it is necessary to write down the equations of motion for the 
set of dynamical variables which in non-Abelian case are x M , p^ and Q a , where Q a is the 
classical colour charge. Unlike in Abelian theory here the colour charge Q a is also subject 
to dynamical evolution if particles interact with external field. The corresponding equations 
of motion are [|37f 



m ^_ = p ^ ( 10 ) 

GST 

m ^ = 9Q a FrPu, (ii) 

dO a 

m^- = -gf abc p»AlQ c . (12) 

Now we can write the collisionless relativistic Boltzmann equation for the one-particle dis- 
tribution fi(x,p, Q) for the i-th type of particle: 



3 Analogous phenomenon is well known for the electron plasma 
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If the external field and the gauge coupling are small the distribution / can be expanded in 
powers of g 

/ = / (0) +^ (1) + -, (14) 

where is the thermal equilibrium distribution function. The induced current which arises 
due to the departure from equilibrium is 

J a = 9 2 Y1 [ dPdQp^i 1 \x,p,Q)Q a} (15) 

where summation stands for all species of particles which are present in the plasma. From 
Eq. dT3p one finds |6J 



= Qa 

and the induced current can be written as 



(16) 



(0) 



J% = 9 2 J2j dPdQp»Q a Q b (A b {x) - ^d Al(x))^-. (17) 



If one takes into account that 



dP = 0- 4 2e(p Q )8(p 2 - 



and 



dQQ Q -|| 5 a6 for quarks (19) 
one arrives at the following expression for the induced current 

J^)=m 2 m j ^vH^d A:( X )-A a (x)), (20) 

where v = (1, p), p = p/|p| and m 2 D0 = l/3g 2 T 2 (N + l/2Nf) for Nf fermion flavours. For 
static fields the above expression simply reduces to 

J*{x) = -m 2 D0 A° a (x)5Z. (21) 

If this expression is inserted to the corresponding field equation the screened Coulomb po- 
tential is easily recovered. 



2.1.3 Polyakov Loop Correlator and the Debye Screening 

Quantities which characterize the phases of a gauge theory are the free energies of static 
configurations of quarks and antiquarks. Following Ref. |]] let us introduce operators, 
ipl(ri,t) and ip a (ri,t), which create and anihilate static quarks with colour a at position r^ 



S 



and time t, along with their conjugates ipl c an d ^ for antiquarks. These static fields satisfy 
equal time anticommutation relations 



[ip a (r i ,t),ipl(r j ,t)] + = SijSab 
The evolution equation for the static quark fields is 

1 d 



^— -t-A q (r i ,t))^(r i ,t)=Q, 



(22) 



(23) 



where f ■ A = t'Mq with r a being the generators of SU(N) gauge group. The above equation 
may be formally integrated to yield 



V(r i5 t) =Texp ^ ^ T-Ao(r h t)j Vfo.O) (24) 
The free energy of a configuration of iV 9 quarks and Nq at positions ri, rjv 8 and r' 1; , r^-_ 



exp(-/3F(ri, ^jy,,^, ,r' 7 



iv„- 



< s|e -/Mf |s >, 



(25) 



with summation over all states with heavy quarks at r 1; rjy, and antiquarks at r[, , r' N _, 

(3 = T _1 and N being the dimension of the SU(N) algebra. The above expression could be 
rewriten as 

exp(-(3F Nq)N -) = 



N N i 



(ab) 



(26) 



where now the sum is over all states with no heavy quarks. Since e ^ generates Euclidian 
time translations, i.e. e^ H 0(t)e~ ,3H = 0(t + (3) for any operator, (KBI) becomes 



exp(-/3F J 



L^E < «'lZ)e-^.i(ri,/3)Vi 1 (ri,0) /?)V>L 9 K, 0) 



(aft) 



x^ 6l (r'^X^, 0) ^ (r^, /3)< 9 _ (r^, 0) |s' > 

Using fl24] ) and ( p2|) and introducing the Polyakov loop 



(f27j) becomes 



L(r) = -TrPexp ( i / dtr • A (r,t) 



exp(-pF Ng>Nq ) 



Tr 



e-^L(r 1 )....LK g )Lt(r' 1 )....Lt(r' / 



A',7 



(27) 



(28) 



(29) 



J This object is also referred to as Wilson line. 
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Figure 1: Two gluon contribution to the Polyakov loop correlation function 

Let us consider the simplest case of an quark-antiquark pair at distance r from each other. 
In this case the aditional free energy due to presence of this pair is simply given by the 
Polyakov loop correlator 

exp(-pF qq )=<L(r)lS(p)>, (30) 

where < .. > denotes thermal average over the gluon configurations. The Polyakov loop has 
following nice properties. First, it is gauge invariant because the gauge field is periodic in 
time 74 M (r,0) = A^r, 0), second, its definition on lattice is straightforward 

1 Nt 

L(r) = -Trl[U (r,t), (31) 
i=i 

where N t is the extension of the lattice in temporal dierction. In the confined phase (T < T c 
with T c being the temperature of the deconfinement phase transition) the large distance 
behavior of the Polyakov loop correlator is the following || 

<L(r)Lt(0)>=exp(- ( r(T)|r|/T) (32) 

implying confinement. In the deconfined phase at very high temperature (small coupling) 
the dominant contribution to the Polyakov loop correlators is given by the exchange of two 
screened electric gluons (see Figure 1). Calculation of the corresponding diagram yields 

<L(r)Lt(0)>=exp(-V(|r|)/T), (33) 

with 

AT 2 _ I / n 2 \ 2 

v u = -nfi-{k) e ~ 2mD °'- (34) 

Thus the screening mass determined from the Polyakov loop correlators is twice the Debye 
mass, Md = 2m p . From definition of V(r) = F q g it is obvious that it is a thermal average 
over qq potentails in the possible colour channels. Let us consider the simplest case of 577(2) 
gauge group. In this case the free energy of quark-antiquark pair (which is often referred to 
as colour averaged potential) can be written as 

e -mr) = 1 (e -WM + 3e -WW) (35) 
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Figure 2: Divergent /-loop contribution to the 2-point function 



where Vi(T4) denotes the qq potential in the singlet (triplet) channel for two isospin | 
particles. These potentials can be calculated perturbatively and contrary to V(r) they are 
dominated by one gluon exchange. The corresponding calculation yields || 

Vi{r) = Ct g 2 ^e- m °o r , (36) 

with C\ = —3/4 and C3 = 1/4. At high temperature g 2 is small and one may expand the 
exponentials in (|35| ) and verify that the 0(g 2 ) terms cancel between singlet and triplet terms 
yielding Eq. (|34"D for V. 



2.2 Infrared divergencies Associated with Static non-Abelian Mag- 
netic Fields 

The application of the naive perturbation theory at finite temperature is obstructed by 
infrared divergencies of thermal field theory. In scalar field theory and in abelian gauge 
theory these IR problems can be cured by appropriate resummations which take into account 
the presence of the screening scale gT (see [p], and also discussion in section 2.5). In non- 
Abelian gauge theories, however, there exists also another sort of IR problems associated 



with static magnetic fields [[17]]. Static magnetic fields are not screened at leading order of 
perturbation theory: the IR limit of the 1-loop transverse self-energy vanishes in all gauges. 

Let us consider a generic /-loop contribution to the self-energy which is shown in Figure 2. 
To evaluate the leading infrared divergent contribution we can assume that all propagators 
in Figure 2 are static, i.e. consider only contribution from zero Matsubara modes. In this 
case neglecting all tensorial structures the contribution of these diagram is 

& (J d3p ) & (p 2 + mT J ^\ (37) 

vertex s v y vertex propagators 
loop integral 

In the above expression we have introduced a magnetic mass since one might expect that 
such a mass scale is generated by higher loop contribution. Using simple power counting 
arguments one can see that for / = 2 the corresponding contribution is 

2 4 T 2 ln- (38) 
m 

and for / > 3 one has 

w^n' 2 (39) 



m 
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Since no magnetic mass is generated at 1-loop level it is natural to assume m ~ g 2 T. In 
this case as one can see from Eq. |3^, infinite number of diagrams contribute to g 4 order 
self-energy. If m is proportional to the power of g, higher than 2 no weak coupling expansion 
could be obtained. However, as we will see in the next section the magnetic mass turns out 
to be of order g 2 T. Different approaches for determination of the numerical value of the 
magnetic mass will be discussed in sections 3 and 4. 

An analysis of the higher loop diagrams similar to that presented above shows that 0(g 6 ) 
contribution to the weak coupling expansion of the free energy is non-perturbative. Let us 
emphasize that the Linde's conjecture does not imply that no weak coupling expansion is 
possible in finite temperature non-abelian gauge theory. It just states that beyond some 
order the coefficients of the weak coupling expansion cannot be evaluated by calculating 
finite numbers of Feynman diagrams. In fact the coefficient of the term proportional to the 
weak coupling expansion of the free energy was calculated by lattice Monte-Carlo technique 
||38|| . Linde's conjecture also implies as we will see in subsection 2.3 that the next-to-leading 
contribution to the Debye mass is non-perturbative. 



2.3 Dimensional reduction 

Finite temparature field theories exist in a volume with one compact dimension, the imag- 
inary time, and the appropriate boundary conditions fl], ||. The extension of the compact 
dimension is j3 = T -1 . At high temparature this dimension becomes arbitrary small. There- 
fore the degrees of freedom in time direction are essentially frozen. This situation can be 
also viewed differently. Consider the Yang-Mills action at finite temperature 



/3 







s = I tit I <i*x\f« v f; v (jo) 



The 4-dimensional gauge field A^(t,x) can be written in terms of an infinite set of 3- 
dimensional fields using Matsubara decomposition 



A£(t,x) = Vf J2 e W ^ n (x), oo n = 2-nT (41) 

n=— oo 

From this expression it is clear that all modes with n ^ have a mass 2irTn. In T — > oo 
limit these modes become infinitely heavy and according to conventional wisdom based on 
Appelquist-Carrazone theorem these will decouple J3§|. Thus the Yang-Mils action reduces 
to an action containing only static fields of the following form 

S = j d 3 x(l^ + l(A^) 2 ), (42) 

with 

D i AZ = d i AZ-g 3 f ab °A$ i A$ (43) 

F% = did} - d 3 A\ + g^ ahc A\ A) (44) 

9l = 9 2 T. (45) 
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This is the action of the 3 dimensional gauge theory which interacts with an Euclidian scalar 
and isovector field A . This matter field originates from the time component of the original 
Yang-Mills field. In more carefull analysis where non-static modes are integrated out at 
1-loop level it has been shown that a mass term and quartic self-interaction is generated for 



the Aq field [53]. Thus the effective action describing the dynamics of zero modes assumes 
the form 



where 



S = J d 3 x(±F«F« + ^(A.Aq) 2 + \m 2 m AlA% + \ A {A%Xtf) , (46) 
m 2 m = l -g 2 T\N+ l -N f ) (47) 



is the leading order Debye mass and 

^A = (G + N-N f )^- 2 (48 
with Nf being the number of fermion flavours [^4], j34| . In Ref. |34| integration of non 



static modes was performed at 2-loop level and it turned out that the 2-loop corrections to 
the parameters g 2 , moo and Xa are small. Higher dimensional and non-local operators are 
also generated by 1-loop integration of non-static modes, however, we have neglected this 
operators because their coefficients vanish in T — > oo limit. The impact of higher dimensional 
operators was analyzed in Ref. JlO| and non-local operators were studied in Ref. ||40|| . Let 



us note that naive decoupling of non-static modes does not take place because the effective 
mass and quartic vertex generated by integration of non-static modes does not vanish in 
T — ► oo limit. If the temperature is sufficiently high, g\ -C m^o i- e - the A field becomes 
very heavy and can be also integrated out. The effective theory obtained by integrating out 
the Aq field is the 3 dimensional Yang-Mills theory 

s YM3 = J (49) 

This theory is superrenormalizable and therefore the only dimensionfull scale in the theory 
is g\. This is the reason why the magnetic mass is expected to be proportional to g 2 = 
g 2 T. Both the 3 dimensional Yang-Mills and the adjoint SU(N) Higgs models are confining 
theories, thus we can see that the infrared problems of finite temperature non-abelian gauge 
theory are related to the physics of 3d confinement. 



2.4 The Non-Abelian Debye Screening beyond Leading Order 

2.4.1 Non-Abelian Debye Screening in One-Loop Resummed Perturbation The- 
ory and the Magnetic Mass 

The next-to-leading correction to IT^ is of order g 3 and arises from summation of ring di- 
agrams The summation of ring diagrams simply amounts to replacing the bare prop- 
agator with the 1-loop resummed propagator and evaluating the 1-loop diagram with the 
resummed propagator. The natural candidate for next-to-leading Debye mass would be 
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ni,(&o — 0, k — > 0). However, this quantity is gauge dependent We have seen in sec- 
tion 2.1 that the Debye screening arises from poles k = ±imr>o in the gluon propagator. 
Therefore the natural extension of the Debye screening mass is 



D 



} k 2 



-m 



Dj 



(50) 



This definition will lead to gauge invariant result because the pole of the propagator was 
proven to be gauge independent at each order of perturbation theory |4"2|| . HL,(k = 0,k) 
was calculated in 1-loop resummed perturbation theory in the high temperature limit in 
Ref. |^3| and it turns out that k 2 = —m 2 D is a starting point of branch cut rather than a 
pole. This fact would mean a non-exponential anomalous screening. If one introduces the 
magnetic mass rriT, i.e. replaces the static massless magnetic propagators with massive ones 



the next-to-leading result for H L in gauge reads 



n L (fc ,fc) 



g 2 TN 
An 



2m 2 D - 2k 2 



rrin 



mo — rn T 
k 



arctan ■ 



(k 2 + m 2 D ) 



k 2 + m 



D 



m\k 



arctan 



m D + rriT 
k 



+ 



m D + V£ m T 



— arctan • 



k 



ttid + tut 



+ 



-4 



(51) 



As one can see from this expression 11^ has gauge invariant poles at k — ±imr) and the 
starting point of the branch cut has been shifted to k = i{mn + Thus the exponential 
nature of the screening is restored. As far as one is interested in obtaining the weak coupling 



expansion of the Debye mass, m D should be replaced by m D0 on the rhs. of equation (|50D . 
Then the next-to-leading correction for the Debye mass reads 



m 



Dl 



— m 



DO 



Ng 2 T 
2vr 



-m D0 



'I - 



rrin 



Am 2 DQ 



In 



2m D0 + m T 
rriT 



1 
2 



rax 
2m D0 



(52) 



Thus we can see that the value of the Debye mass at next-to-leading order is strongly 
influenced by the magnetic mass. An alternative possibility which will be considered in 
section 3 is a self-consistent determination of from the equation flollD, which is equivalent 
diagrammatically to a summation of "super-daisy" diagrams [f4"4 |. Analogously to the Debye 
mass the magnetic screening mass can be defined as 



ran 



rir(fco = 0, k 2 = —m 



T)i 



(53) 



where I1 T is defined by equation (||). Though the pole of the transverse gluon propaga- 
tor is called magnetic screening mass its relation to the screening of static magnetic fields 
is not obvious. The magnetic screening could be defined in terms of the free energy of 
monopole-antimonopole pair 0. However, it cannot be stated that for the weak coupling 

The divergence of the static non-Abelian magnetic field is non-zero because the the corresponding non- 
Abelian field equation is V • B a = gf abc A b ■ B c . Therefore there may exist non-Abelian field configurations 
which act as magnetic monopoles 
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Figure 3: Contribution of magnetic gluons to the Polyakov loop correlation function 



regime this energy is dominated by exchange of a static transverse gluon. Because of the 
Linde's conjecture all multi gluon exchange have to be considered. 

The magnetic screening can be defined by studying the large distance behaviour of the 
field created by a magnetic monopole. This was realized on lattice investigations [45, [46 



by calculating the magnetic flux leaving the spatial 3d cube and also by semiclassical ap- 
proximation by Bird and Miiller |f4"7| . In the lattice investigation the corresponding inverse 
screening length was found to be 0.24g 2 T [4~5| and 0.27(3)g 2 T in Ref. P6] |. The semiclas- 
sical approach of Ref. E_7| yields the result 0.25g 2 T which is consistent with the previous 
ones. From the arguments given above it is clear, however, that this inverse screening length 
cannot be directly related to the magnetic mass defined from the propagator. 



2.4.2 Non-Perturbative definitions of the Debye mass 

We are interested in definitions of the Debye mass which do not rely on perturbation theory 
and can be easily implemented in lattice Monte-Carlo simulations. A natural way to deter- 
mine the Debye screening mass non-perturbatively is to measure the gluon propagator on 
lattice. Such measurement were done in Refs. ESI, HU[ in finite temperature 4d SU(2) gauge 



theory and also in 3d effective theory in Ref. |50| . We will discuss the non-perturbative 
determination of the propagators in section 4 in detail. 

Let us consider the Polyakov loop correlator which is a manifestly gauge invariant quan- 
tity. At leading order the Polyakov loop correlator is determined by the exchange of two 
electrostatic gluons as shown in Fugure 2 (see subsection 2.1.3) leading to exponential screen- 
ing with the inverse screening length equal to 2rriDo- 

Beyond leading order, however, Polyakov loops can also couple to magnetostatic gluons 
as shown in Figure 3. At high temparatures different mass scales which are present in the 
plasma are separated 2irT ^> gT ^> g 2 T due to assymptotic freedom. Then the infrared 
problems associated with static magnetic fields are related to the physics of 3 dimensional 
confinement. Because of the 3 dimensional confinement Polyakov loops cannot exchange a 
pair of magnetic gluons, the pair will instead form a glueball with a mass of order g\ = g 2 T. 
Thus the large distance behaviour of the Polyakov loop correlators will be determined by the 



lightest glueball state of the 3 dimensional Yang-Mills theory [[51], Thus the true large 
distance behaviour of Polyakov loop correlator at very high temperature has nothing to do 
with the physics of the Debye screening. 

For vectorially coupled theories (like QCD at zero chemical potential but not the elec- 



15 



troweak theory) an alternative non-perturbative definition of the Debye mass was suggested 



by Arnold and Yaffe This definition relies on the fact that in vectorially-coupled theories 
at finite temperature there is a symmetry, namely the Euclidian time reflection symmetry 
that excludes the unwanted exchange of pairs of magnetostatic gluons. The crucial property 
is that the temporal component of the gauge field Aq is odd under this symmetry while the 
spatial gauge field Aj is even. Thus non-perturbatively the Debye mass can be associated 
with a correlation function of some gauge invariant operator which is odd under Euclidian 
time reflection and which has the largest correlation length. Some gauge invariant time- 
reflection odd operators are listed in Ref. P5| . As far as one is interested in the leading 
non-perturbative correction i.e. the coefficient of the g 2 lng and g 2 terms in 

m D = m D0 + dg 2 \ng + c 2 g 2 + ... (54) 

it is not necessary to calculate correlation function in the full 4 dimensional theory. One can 
instead consider the effective theory, construct the corresponding operators and measure their 
correlation function there. Various realisation of this strategy have been recently discussed 
in the literature 13(1 |3Tj, |32], |33 . 



2.5 Screened Perturbation Theory 

One of the main motivation for studying the screening masses is their role in the physical 
expression of different thermodynamic quantities like thermodynamic potentials. The weak 
coupling expansion of the QCD free energy does show very bad convergence properties 

Only in the TeV temperature range can one find a satisfactory 



M, 55 



57, 58 



numerical convergence rate . 
Recently, in the case of scalar 



!> 4 theory, where conventional perturbative expansion 
encounters similar convergence problems, several proposals have been made for new types 



of perturbative expansion which shows better convergence |HT|, . The outcome of the 
resummed perturbative calculation for the free energy density may depend on the correct 
choice of the screening mass pO, 
In Ref. 



60 



it was shown that the poor convergence of perturbative expansion can be 
improved if the series expansion is reorganized in loop expansion with screened propagators; 
i.e. instead of expanding around the free energy density of a massless ideal gas, loop expan- 
sion is performed starting from a massive ideal gas |. Following Ref. |60| consider an O(N) 
symmetric scalar field theory with the following Lagrangian 



L 



L — Lq + Li n t, 

= \m 



2 + m 2 d> 2 ) 



L 



int 



im 2 0?+2w(^r + ^(^-l)(^ 2 



2 

247V 1 



&\ 2 



(55) 



Following ref. |)3], we have introduced and substracted a mass term with m ~ O(g) 
in order to reorganize the perturbative expansion. The coupling constant renormalisation 



factor is given by |65 



3g 2 N + 8 
(An) 2 18e 



0(g 4 



(56) 



6 The propagators become massive by the mechanism of Debye screening generated by the thermal fluc- 
tuation in analogous way as in gauge theories. 
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Figure 4: Shown is the free-energy density of the N=l scalar field theory as a function of 
the scalar self-coupling g in the units of the free energy density of a massless ideal gas ( ZL ^j-) 
(from Ref. J60|). The curves represent 2- loop (a) and 1-loop (b) results of the screened loop 
expansion as well as the 0(g 2 ) (c) and 0(g 3 ) (d) results of the conventional perturbative 
expansion. 



and we have neglected the field renormalisation factor Z%, which is unity up to 0(g 4 ). 

The free energy can be easily evaluated with massive propagators up to 2-loop, but in the 
limit of large N also the 3-loop contribution can be evaluated. The details of the calculation 
can be found in Ref. |6(J . To evaluate the free energy numerically the screening mass has to 
be specified. It turns out that the final result is sensitive to the actual value of the screening 
mass [|66| . The most natural way of the determination of the screening mass is the use of a 
1-loop gap equation [[44], [67], |68| , which for our case reads [] 



m 



g 2 (N 



6N 



Mm) + 



m 



16tt 2 



with 



1 171 ) 



d 3 p 



(27T 



i 3 p 2 + u 2 + m 2 ' 



2nTn. 



(57) 



(5* 



The numerical results for the free energy density in the case of the N — 1 component 
theory as a function of the self-coupling calculated in the screened perturbation theory is 
shown in Figure 4. Also shown are the results obtained with the conventional perturbative 
expansion up to 0(g 2 ) and 0(g 3 ) in the scalar self-coupling. The Figure clearly shows that 
the difference between the 1-loop and 2-loop results of the screened perturbative expansion 
is small and essentially stays constant over the entire range of couplings. In the conventional 
perturbative approach, on the other hand, different orders in the expansion differ largely 
from each other. For large N as it was mentioned above also the 3-loop contribution can 
be calculated. The corresponding results are shown in Figure 5 for iV = 10. As one can see 
from the figure the convergence of the perturbative expansion is substantionally improved 
as compared to the usual perturbative expansion. 

Let us notice that an attempt to generalize this approach to the case of QCD have been 
recently reported in Ref. |69| The authors were able to find a compact formula for the 



7 In the limit of large iV this equation is exact and therefore its solution can also be used in 3-loop 
calculations. The gap equation should be properly renormalized before it is solved numerically. Therefore a 
divergent and scale dependent counterterm was added; see discussion in Ref. [p2] on this issue. 
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Figure 5: The free energy density of the scalar field theory in the large- N limit as a function 
of the scalar self-coupling g in units of the free energy density of a massless ideal gas (from 
Ref |60]]). Shown are results for N = 10 (see text). The curves represent 3- loop (a) 2-loop 
(b) and 1-loop (c) results of the screened loop expansion as well as the 0(g 4 ) (d) results of 
the conventional perturbative expansion. 

leading order massive gas contribution to the free energy density of the hot gluon gas. In 
distinctions to phenomenological approaches |T5|] , their expression also takes into account 
the contribution of the longitudinal fluctuations and also the effect of Landau-damping 
appears in addition to the quasiparticle piece of the free energy. The screening mass of 
the gluons remains a free parameter, which emphasizes the importance of its independent 
determination. 
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3 Coupled Gap Equation for the Screening Masses in 
Hot SU(N) Gauge Theory 

As it was discussed in section 2.4 a natural way to calculate the Debye mass is to solve the ( 
gap ) equation m 2 D = n i (fc = 0, k 2 = —m 2 D ) in a self-consistent way This definition is gauge 
invariant but requires the introduction of the magnetic mass which can be associated with a 
dynamically generated mass gap of pure 3d gauge theory if the temperature is high enough. 
The magnetic mass was calculated using self-consitent approach in the effective 3d gauged 



Higgs and also in 3d pure gauge theories q7|, [7(], [72|, [731 . In the case of the Electroweak 
Theory (or SU(2) Higgs model ) the gauge coupling is small for the temperature range of 
interest and a hiearchy of the different mass scales holds 2ttT » gT » g 2 T. Due to this 
fact different screening masses could be determined separately in the corresponding effective 
theories. In hot SU(N) theory for the temperature range of interest the above hierachy of 
scales fails to hold. Therefore we should investigate a coupled set of gap equations for all 
the screened modes and determine the corresponding screening lengths simultanously. 

The most straightforward way to do this would be to derive gap equations in the full four- 
dimensional theory. However, it is not known how to generalize the by-now well-established 
three-dimensional gauge invariant resummation J70, [H], |72], |n| to four dimensions. There 



are gauge transformations which might mix static and non-static modes, therefore the re- 



summation of the static modes only, which was suggested in [|74| violates gauge invariance. 
Since screening masses are static quantities it is natural to calculate them in the framework 
of an effective three-dimensional theory which is, however, valid only up to scales k < gT. It 
was shown that the accuracy of the description of such theories is improved if in the action of 
the effective theory besides superrenormalizable operators one also takes into account higher 
dimensional and non-local operators |7f| [40[| . 

When deriving the gap equation in the framework of three-dimensional adjoint Higgs 
model it is important to address the question whether the symmetric or the broken phase is 
the physical one. Dimensional reduction is valid if A -C y/T. However, in the broken phase 
the expectation value of the Aq field is Aq ~ ~. Therefore for small g the broken symmetry 
phase cannot be the physical one. We will return to this question in section 4, but here we 
will assume no A background. 

In this section we shall assume only the separation of the static and non-static scales 
in the pure SU(N) gauge theory, that is we assume that g(T) « 2tt, and derive different 
versions of the coupled gap equations in the corresponding effective theory, the 3d adjoint 
Higgs model. Presentation in this section relies on Ref. 



3.1 Gauge non- invariant resummation scheme 



This scheme for the evaluation of the magnetic mass was first suggested in WE and for the 



Debye mass in |43| . It was noticed in Ref. JFlJ] that magnetic mass obtained in this scheme 



is gauge dependent and therefore cannot be regarded as physically meaningfull. However, 
even in gauge invariant resummation schemes the value of the magnetic mass cannot be 
defined unambigously, because its approxiomate value depends on the specific resummation 



scheme [73]. Therefore it is of a certain interest to calculate the magnetic mass even in 



a gauge non-invariant scheme just to compare the amount of gauge dependence with the 
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ambiguity of different gauge invariant approaches, when evaluated to leading order is some 
approxiomation scheme. 

The three-dimensional effective Lagrangian relevant for 1-loop calculations can be written 
with i?g gauge as (see Eq. (£f6|)) 

L = \F%F§ + \ ((A^o) 2 + mlA-Al) + 
\m\AlAl + (diC a D iC a + m 2 G c a c a ) + ^(d^) 2 + L ctl 
L ct = \{m 2 m - m 2 D )AlAl - \ni%.A$A$ - m 2 G c a c a (59) 

Here we have added and substracted a mass term for A iy A and the ghost fields (c) with 
their exact screening masses. m 2 D0 is the tree-level (from the point of view of the effective 
theory) Debye mass for Aq, which was generated during the procedure of the dimensional 
reduction and defined by Eq. (f47|). The propagators can be read from the quadratic part of 
the Lagrangian and can be found in the Appendix A, where one also finds some details of the 
evaluation of the relevant Feynman diagrams contributing to the different 2-point functions. 
It should be noticed when performing the resummation of the pure gauge sector with the 
unique mass term my, the longitudinal and transverse gluons acquire different masses, which 
are, however related by = y/^rriT {ttit is the transverse and is the longitudinal mass). 
It is also possible to perform the resummation by introducing independent masses for the 
longitudinal and transverse gluons, but then the corresponding gap equations will have only 
complex solutions. The gauge boson self-energy can be decomposed as 

Ilij(k) = (Sij ^)II T (fc, m T , m D , m G ) + -j^-Tl L (k, m T , m D , m G ). (60) 

In the Appendix A we give the expression of n^- in terms of a few fundamental three- 
dimensional loop-integrals. These integrals are easily evaluated and an explicit but very 
cumbersome functional form can be written for the longitudinal and transversal projections 
of the polarisation matrix. The self-energy of the ghost fields is also given in Appendix A. 



The self energy for Aq was first calculated in [|43f| : 



n o(A;, m D , m T ) = m 2 m + 



2 , g 2 N 



Lit 



, 2(m 2 n -k 2 -m% 12) 

-mo — fn r p H — — r — L1 - L arctan ■ 



k rriD+rriT 



+ (k 2 + m 2 D ) ^(arctan - arctan ^) + (^ - 1)^ j . (61) 



The above self-energies are functions of the magnitude of momentum but also depend 
parametrically on the screening masses. The on-shell gap equations now can be written as 

m\ = Hr(fc 2 = — m^, mr, JTic, wg), 

m 2 D = II 00 (A; 2 = -m 2 D ,m D ,m T ), (62) 
m 2 G = Il G (fc 2 = -m G ,m T ,m G ). 

On the mass-shell n o(& 2 = —ml^mo^rriT) is gauge parameter independent, but Xl T {k 2 = 
—m 2 ^,mT,'mD,'m G ) and H G (k 2 = —m G ,mT,m G ) do depend on the gauge fixing parameter, 
therefore the masses obtained from this coupled set of gap equations are gauge dependent. 



20 




Figure 6: The temperature dependence of m D /m D0 (solid line) and m T /m pp for Feynman 
(£ = 1) gauge, moo is the leading order result for the Debye mass and m pp is the value of 
the magnetic mass obtained by Buchmuller and Philipsen for pure SU (2) gauge theory. 
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Figure 7: The dependence of mp /moo (solid line) and rriT/m pp on gauge parameter £ at 
T = 10 4 T C , m D0 is the leading order result for the Debye mass and m pp is the value of the 
magnetic mass obtained by Buchmuller and Philipsen for pure SU (2) gauge theory. 
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In the following numerical investigations we shall consider the case of the SU(2) gauge 
group. The 4- dimensional coupling constant is taken at scale \x = 2ttT, where /x is the MS 
scale and 1-loop relation for the gauge parameter of the effective theory is used. To set the 
temperature scale we use the relation T c / Kjjg = 1.06 obtained from numerical simulation of 
the finite temperature SU(2) gauge theory |4"9|| . 

The temperature dependence of mo in the Feynman gauge is plotted in Figure 6. As 
one can see from the plot receives 30% positive correction compared to the leading 
order result, while the magnetic mass stays very close to the value calculated by Buchmiiller 
and Philipsen [[7TJ, given below in Eq.(^). Since the masses are gauge dependent in this 
approach, it is important to investigate the dependence of the screening masses on the gauge 
parameter £. It turns out that one gets real values for the masses from the gap equations 
only if £ 6 [1>5). The dependence of the screening masses on £ in this range is shown in 
Figure 7. One can see, that the ^-dependence of tut in this interval is 40 %, while for mo 
it remains in the 10 % range. 



3.2 Gauge Invariant Approach 

Gauge invariant approaches for the magnetic mass generation in three-dimensional pure 
SU(N) gauge theory were suggested by Buchmiiller and Philipsen (BP) [[F(J and by Alexanian 
and Nair (AN) f72|. The basic idea of both approaches is to add and substract a gauge 
invariant mass generating term to the original Lagragian: 

L = lp*F% + 5L - 15L, (63) 

where / is the (formal) loop expansion parameter. In the final result one should set I = 1 of 
course. In the scheme of BP SL is just the Lagrangian of the gauged non-linear sigma model. 
In the scheme of AN the mass generating term relevant for 1-loop calculation assumes the 
form [72 1 

5L = ~ ^)A j + f abc V l]k A1A)A^ (64) 



where the explicit expression for can be find in Ref . |72| . Till now only these two gauge 
invariant schemes are known to provide real values for the magnetic mass ||73|| . 

In these approaches the gauge boson self-energy is automatically transverse and there 
is no need to project the transverse part from the polarisation tensor. The corresponding 
expression for the on-shell self-energy reads 

Uxik = imT,rriT) = CrriT, (65) 

where 

IS! >n 3 -I], BP. <66) 

In gauge invariant approaches the ghost field remains massless and there is no gap equation 
corresponding to it. 

Since we are interested in calculating the screening masses in the three-dimensional SU(N) 
adjoint Higgs model, IIt(A;, tut) should be supplemented by the corresponding contribution 
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coming from A fields. This contribution is calculated from diagrams d) and e) of the 
Appendix to be 

XT1 A on \ 9 2 N ( m D k 2 + Am 2 D k \ f k { kj\ 
5Il^ (k, m D ) = — — — — H — - arctan -j- . (67) 



v v ' An V 2 4k 2m D J \ J k 2 

It is transverse and gauge parameter independent, it also does not depend on the specific 
resummation scheme applied to the magnetostatic sector. It should be also noticed that it 
starts to contribute to the gap equation at 0(g 5 ) level in the weak coupling regime, thus 
preserving the magnetic mass scale to be of order g 2 T. This is the reason why no "hierarchy" 
problem arises in this case, at least for moderate g values. 

The self energy of A depends on the specific resummation scheme. To calculate this 
one has to fix a specific gauge. In Ref. J73[ it was shown that it is possible to integrate 
out the sigma fields and the resulting Lagrangian describes massive gauge fields and has no 
gauge-fixing and ghost terms. Therefore to calculate the self-energy of Aq in BP scheme one 
should evaluate diagrams / and g in appendix A with massive vector boson propagators. 
The corresponding calculation leads to 



n 00 (fc, m D , m T ) = m 2 m + 



9 2 N 
4tt 



, 2<m 2 n -k 2 -m 2 /2) 

-mo — tut H — - 12 — r — L1 - L arctan 



rriT+rno 



(/2+ ^(-m T + ^ arctan^)]. (68) 



This expression is different from the expression of n o calculated in the gauge non-invariant 
approach (see eq. (|61| ) ) but its analytic properties and on-shell value is the same as of (|6"ID . 
For the resummation scheme of AN it is important to find a convenient gauge fixing term 
which reads 

L ff=\ d ^{ l - m T^)d 3 A« (69) 

With this gauge fixing the expression for the self-energy of A coincides with (|6T| ) if it is 
evaluated at £ = 1. The coupled set of gap equations now can be written as 

= trriT, Trip), 

m 2 D = U 00 {k = im D ,m D ,m T ). (70) 

The temperature dependence of m D obtained from this coupled set of gap equations is shown 
in Figure 8 for both schemes, where we have again normalized the Debye mass by the leading 
order result, moo- The temperature dependence of the magnetic mass is shown in Figure 9, 
where we have normalized m^ by the value of the magnetic mass obtained for pure three- 



dimensional SU(2) theory, in the BP (AN) gauge invariant calculations 0, 72]. As one can 
see the contribution of A to the magnetic mass is between 1 and 3%. From Figures 8 and 
9 it is also seen that the temperature dependence of the screening masses is very similar to 
the temperature dependence of the respective leading order results. 



3.3 Contribution of non-Local Operators to the Gap 
Equation 

In the previous section the gap equations were derived for an effective local superrenor- 
malizable theory. In this case the effect of non-static modes in the 2-point function were 
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Figure 8: The temperature dependence of the scaled Debye mass for BP resummation scheme 
(solid) and for the AN resummation scheme (dashed). The scaling factor is moo- 




Figure 9: The temperature dependence of the scaled magnetic mass for BP resummation 
scheme (solid) and for the AN resummation scheme (dashed). The scaling factors are m^ p 
and m^ N , respectively. 
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represented by the thermal mass for A and by the field renormalization factors which relate 



3d fields to the corresponding 4d ones. It was shown in |TT, 77]] that when one performs the 
procedure of the dimensional reduction in gauges the parameters of the superrenormaliz- 
able effective theory are gauge independent, only the expressions of 3d fields in terms of 4d 
ones depend on the gauge parameter. However, this does not hold for higher dimensional 
and non-local operators, which are generally gauge dependent. At 1-loop level the only 
diagrams contributing non-locally to the gap equation are those which have two non-static 
line inside the loop, diagrams with one static and one non-static line inside the loop are 
forbiden because of 4-momentum conservation. Therefore at 1-loop level the non-locality 
scale (27rT) _1 is much smaller than the relevant length scales. 

In general the non-static contribution to the static 2-point function U^ u (ko = 0,k) can 
be written as 

An-^o = o, k) = syMk) + slHfc - ^)n(fc), 

ir»(*0 = &!(/*, Ok 2 + r 2 E: =2 &n(0(A) 2n , (7i) 

where the coefficents a n and b n can be calculated for arbitrary n. The first two terms in the 
expressions of Il 00 and first term of II(fc) are already included into the 3d effective theory as 
part of the tree level mass and the definition of the 3d fields in term of 4d fields. The last two 
sums will contribute to the 3d effective lagrangian as quadratic non-local operators. There 
are also higher dimensional operators as well as non-local 3- and 4-point verticies in the 
effective lagrangian, however, since we restrict our interest to 1-loop gap equations these are 
not important for us. Their contribution would correspond to 2 or higher loop contribution 
in the full four-dimensional theory. 

We have estimated by direct numerical evaluation of the infinite sums in ([71]) the contri- 
bution of non-local operators to be less than 1% in the temperature range T = (3 — 10 4 ) T c , 
thus neither their contribution nor their gauge dependence is essential. 



3.4 What Have We Learnt from the Coupled Gap Equations ? 

In this section an attempt have been made to extract the electric and magnetic screening 
masses from the coupled set of gap equations of the three-dimensional SU(N) adjoint Higgs 
model considered as an effective theory of QCD. The screening masses have been studied 
using gauge non-invariant as well as gauge invariant resummation schemes. In the gauge non- 
invariant formalism we have observed rather strong gauge parameter dependence, therefore 
the results extracted from it are not very informative. It is still interesting to note that in 
Feynman gauge (£ = 1) the results for the magnetic mass are rather close to those obtained 
from the gauge invariant resummation scheme of Buchmuller and Philipsen. In gauge invari- 
ant treatments we have compared two different resummation schemes, that of Buchmuller 
and Philipsen (BP) and one proposed by Alexanian and Nair (AN). Qualitatively these two 
resummations lead to similar results, but in the BP scheme one has smaller magnetic mass 
and larger Debye mass than in the AN scheme. 

Let us summarize our view on the interaction of the electric A, and the magnetic Ai fields. 
In both schemes one can see that the dynamics of Aq is largely influenced by the magnetic 
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sector, however, no similar feedback on the magnetic sector is seen. The magnetic masses 
calculated from the coupled gap equations provide screening masses which are 1% smaller 
than evaluated in the pure gauge theory. This fact suggests that the influence of the adjoint 
scalar field is similar to that of the fundamental Higgs field, because the magnetic mass 
calculated in the symmetric phase of 3d SU(2) Higgs theory by the gap equation technique 
is roughly the same as in pure gauge theory [[7(| . 

Finally we compare our results with recent Monte-Carlo data for the screening masses 
obtained in 4d finite temperature SU(2) gauge theory |]¥9[| . The data on the magnetic mass 
found from this simulation in the temperature range T = (10 — 10 4 )T C can be fitted well by 
the formula my = 0.456(6)g 2 T. This value of the magnetic mass is considerably larger than 
what one obtains from the magnetic gap equation and 3d simulations, where the results are 
approximately mx = 0.28g 2 T for the BP scheme, 0.38g 2 T for the AN scheme and 0.35(l)g 2 T 
from 3d simulations. The corresponding results for the 3d adjoint Higgs model are also close 
to these values. 

The data from Monte-Carlo simulation for the Debye mass in the above mentioned tem- 
perature interval can be fitted using the following leading order-like ansatz a/ 1.69 (2) g (T)T 
|49|| , which means that the Debye mass is roughly 1.6m£>o, where mo is the leading order 



result. The gap equations at the same time, as one can see from Figure 8 give (1.2 — 1.3)mooi 
depending on the resummation scheme. While there is no quantitative agreement between 
masses measured in Monte-Carlo simulation and those obtained from the gap equation, the 
temperature dependence of these masses in the temperature interval T = (10 — 10 4 )T C seems 
to follow the temperature dependence of the leading order result. 

Finally let us discuss the consistency of the gap equation approach. The validity of 1- 
loop gap equation for the magnetic mass in pure gauge theory was critically questioned in 
Refs. [78, [TjJ. We have seen that the two different resummation schemes lead to somewhat 
different results. If the gap equation aproach is valid the difference between the values of 
the magnetic masses calculated in two schemes — m^ p should be of the same order in 
magnitude as the 2-loop correction to the magnetic mass [73|. A recent 2-loop analysis of 
the magnetic mass yields = 0.34gf |80| . Thus the gap equation approach seems to be 
viable. Moreover in Ref. [pl| a strong coupling expansion for the magnetic mass was carried 



out with the result 0.32g§. In next section we will discuss the determination of the magnetic 
mass using lattice Monte-Carlo technique. 
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4 Screening Masses of Hot SU(2) Gauge Theory from 
Monte-Carlo Simulation of the Lattice 3d Adjoint 
Higgs Model 

Since screening masses are static quantities it is expected that they can be determined in 
a 3d effective theory, the 3d SU(N) adjoint Higgs model, provided the temperature is high 
enough. However, in the case of the hot SU(N) gauge theory one may worry whether the 
standard arguments of dimensional reduction apply. First of all the coupling constant is large 
g ~ 1 for the physically interesting temperature range and thus the requirement gT << ttT 
is not really satisfied. Another problem is the identification of the physical phase of the 
3d adjoint Higgs model, i.e. the phase which corresponds to the high temperature phase 
of SU(N) gauge theory 0. The 3d adjoint Higgs model is known to have two phases, the 
symmetric (confinement) phase and the broken (Higgs) phase [[52], |S3], |34J separated by a line 
of strong I s * order transition for small Higgs couplings. The perturbative calculation of the 
effective potential p|, EM suggests that the Higgs phase of the reduced theory corresponds 



to the deconfined phase of the 4d SU(N) theory. This conclusion seems to be supported by 
the 2-loop level dimensional reduction and the analysis of the phase diagram performed in 
Ref. pi| . Though an earlier study [^] indicated that the symmetric phase of the 3d adjoint 



Higgs model is the physical one, a more carefull analysis shows that this conclusion may be 
altered by higher order corrections generated in the process of dimensional reduction (see 
appendix B.l). The fact that the broken phase turns out to be the physical one is contrary 
to naive expectations (see the discussion in the previous section). This circumstance and the 
0(1) value of the gauge coupling raises serios doubts about the validity of the dimensional 
reduction. Therefore a non-perturbative study is nedeed to clarify the applicability of the 
3d effective theory. 

The aim of this section is to clarify whether the screening masses, defined as poles of the 
corresponding lattice propagators can be determined in the effective theory for the simplest 
case of the SU(2) pure gauge theory. For this model precise 4d data on screening masses exist 
for a huge temperature range Inclusion of fermions into the model is straightforward f\ 
and generalization to the SU(3) case will not introduce any novel feature. 

Gauge invariant definitions of the Debye mass as it was discussed in section 2 rely on the 



3d effective theory [^, |30], |34|. Although these definitions are the best in the sense that they 
are explicitly gauge invariant, little is known about the corresponding correlators in the full 
4d theory. Gauge invariant correlators yield larger masses than the pole masses |34| and 



further studies are required to establish the connection between them. We will try to relate 



these masses in section 4.5 of the thesis in the spirit of the so-called additive model [IT3 . 



8 In the vicinity of the deconfinement phase transition the gauge coupling is very large and the dimensional 
reduction is certainly not valid. Therefore the two phases of the 3d adjoint Higgs model has nothing to do 
with the phases of 4d SU(N) gauge theory 

9 Fermionic fields do not appear in the effective theory, their role is only to modify the parameters of the 
effective theory. 
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4.1 The 3d SU(2) Adjoint Higgs Model on Lattice 

The lattice action for the 3d adjoint Higgs model used in the present paper has the form 

S = P E \ TrU P + P E ^o(x)^(x)i (x + i)U}(x) + 
p 



X,l 



E 



l.\ 1 r „ , „ (\ ~ 



-P(3 + -h -Tr^(x) +0x1 -Tr^(x; 



(72) 



where t/p is the plaquette, JJ\ are the usual link variables and the adjoint Higgs field is 
parameterized, as in [[34]] by anti-hermitian matrices A = i (y a A^ (a a are the usual Pauli 
matricies). Furthermore /3 is the lattice gauge coupling, x parameterizes the quartic self 
coupling of the Higgs field and h denotes the bare Higgs mass squared. The details of the 
derivation of the lattice action can be found in appendix B.l. 

In principle the parameters appearing in eq. (|72| ) can be related to the parameters of 



the original 4d theory via the procedure of dimensional reduction [33] which is essentially 
perturbative and discussed in appendix B.l. The obvious problem with this approach, as it 
was already discussed before, is that the physical phase turns out to be the broken phase. 
This problem may signal that higher dimensional and non-local operators have to be kept 
in the effective action to ensure the consistency of the whole approach Certainly it is 
difficult to take into account all possible (higher dimensional and non-local )operators in a 
lattice investigation. 

There are two ways out from this unfortunate situation. One can assume that the com- 
plicated action containing all kinds of higher dimensional and/or non-local operators can be 
mapped onto the effective action (|72|) which is the action of some coarse-grained theory with 
some finite UV cutoff and does not necessarily correspond to the naive continuum 3d adjoint 
Higgs model. Then parameters appearing in ( ff2| ) are found by matching some quantities 
which are equally well calculable both in the full 4d lattice theory and in the effective 3d 
lattice theory, i.e. by a non-perturbative matching. 



Another possibility which was suggested in Ref. [|34j is based on observation that the 



parameters of the 3d adjoint Higgs model which correspond to the high temperature SU(2) 
gauge theory lie in the region of the parameter space which is close to the transition line and 
where for any finite volume available in numerical simulations the symmetric phase turns 
out to be metastable (see appendix B.l for a detailed discussion on this issue). For practical 
purposes, e.g. measurements, the metastable phase behaves as if it was stable and thus 
parameters appearing in (^) can be fixed at their values given by standard dimensional 
reduction. One proceeds futher by arguing that if higher dimensional and/or non-local op- 
erators were taken into account the transition line would be shifted so that the parameters 
corresponding to 4d physics would move to the region where the symmetric phase is abso- 
lutely stable but, the inclusion of these operators would have little impact on the values of 
the quantities we are interested in. 

Here we will consider both possibilities. Comparison between corresponding 3d and 4d 
data, hopefully enables us to clarify which scenario takes place. 

10 An attempt to caclulate the higher dimensional operators in the effective 3d action was recently done in 
Ref. M. 
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Let us start to discuss the first strategy for fixing the parameters j3, h, x, namely the 
non-perturbative matching. We will turn to the discussion of the second strategy in the 
next subsection. The deconfined high temperature phase of the 4d SU(2) gauge theory 
corresponds to some surface in the parameter space (j3, h, x), the surface of 4d physics. For 
fixed value of f3 the 4d physics is described by a line on this surface, which will be refered to 
as the line of 4d physics. One of the aim of numerical simulations is the determination of the 
surface of 4d physics. In general this would require a matching analysis in a 3d parameter 
space (P, x, h), which is clearly a difficult task. We thus will follow a more moderate 
approach and fix two of the three parameters namely (3 and x, to the values obtained from 
the perturbative procedure of dimensional reduction. The values of these parameters at 
2-loop level are |34 



P 
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i > 



- 9 2 (M) T 
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g 2 (v) 
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^ln^^ 

3 7.0555T' 



(73) 
(74) 
(75) 



with a and T denoting the lattice spacing and temperature, respectively. The coupling 
constant of the 4d theory g 2 (/i) is defined through the 2-loop formula 
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(76) 



The parameter h will be left free to allow non-perturbative matching. In order to be able 
to compare the results of the 3d simulation with the corresponding calculations in the 4d 
theory it is necessary to fix the renormalization and the temperature scale. We choose the 
renormalization scale to be jj, — 2ttT, which ensures that corrections to the leading order 
results for the parameters g\ and x of the effective theory are small. Furthermore we use 
the relation T c = 1.06A^gg from |4"9fl . Now the temperature scale is fixed completely and the 



physical temperature may be varied by varying the parameter x. The lattice spacing was 
chosen according to the criterium a « rrT 1 « Na, where N is the extension of the lattice 
and m is the mass we want to measure. 



4.2 Determination of the Propagators and the Screening Masses 
on Lattice 

The main goal of this section is to study the propagators of scalar and vector (gauge) fields. 
For this purpose one has to fix a specific gauge, which is chosen to be the Landau gauge 
d^A^ = 0. The choice of this gauge is motivated by the fact that the 4d Landau gauge 
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condition used in [^5, [RJ for static field configuration is equivalent to the 3d Landau gauge 
condition. On the lattice this condition is realized by maximizing the quantity []36 |: 



TV 



£([/,(x) + [/t( x ; 



X,2 



(77) 



The gauge fixing is performed using the overrelaxation algorithm [87]], which in our case is 
as efficient as the combined overrelaxation and FFT algorithm used in ff^| . The vector 
field is defined in terms of the link variables as 



I,(x) = i(Oi(x) - ltf(x)) 



(78) 



We are interested in extracting the electric (Debye) mp and the magnetic uit screening 
masses from the long distance behaviour of the scalar and vector propagators defined as 



G D (z) =< TrA Q (z)Al(0) >~ exp(-m D z), 
Gt{z) = -{G x {z) + G 2 (z)) ~ exp(-m T z), 



(79) 
(80) 



with 



Gi(z) =< Tr Ai(z)Ai(0) >, 

^A^(x,y,z), /i = 0, 1, 2 



Note that due to the Landau gauge condition G^(z) should be constant. This fact can 
be used to test the precision and validity of the gauge fixing procedure. In our case this 
condition is satisfied with an accuracy of 0.001%. 

Besides the scalar and vector propagators of the adjoint Higgs model we also calculate 
the gauge invariant scalar correlators and analyze the propagators also in the limit of a 3d 
pure gauge theory. 

To extract the masses from the correlation functions we have used the general fitting 
ansatz 



A 



exp(-mz) exp (— m (N z — z)) 



+ 



(82) 



Previous investigations of gauge boson and quark propagators in Landau gauge PBL R8L B9| 



have shown that effective masses extracted from the correlation functions rise with increasing 
Euclidean time separation and eventually reach a plateau. This is contrary to local masses 
extracted from gauge invariant correlation functions which approach a plateau from above 
and is a direct consequence of a non-positive definite transfer matrix in Landau gauge. A fit 
with b 7^ allows to extract stable masses already at shorter distances. Most of our numerical 
studies have been performed on lattices of size 32 2 x 64 0. From previous studies of gluon 

11 Additional calculations also have been performed on a 16 2 x 32 lattice and it have been checked explicitly 
that results for the electric mass show no volume dependence. Furthermore some calculations have been 
performed on a 32 2 x 96 lattice to check that the lattice size used for our calculations was sufficient for the 
determination of the magnetic mass. 
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propagators in the 4- dimensional SU(2) gauge theory we know that such large lattices are 
needed to observe a plateau in local masses [(49)]. In particular, this is the case for the rather 



small magnetic screening mass which leads to a rather slow decay of the correlation functions. 
We thus use the same spatial lattice size as in those studies. 

We have used correlated (Michael-McKerrell) fits with eigenvalue smearing |Mj . Our fits 
have been constrained to the region where local masses show a plateau (typically z ~ 15 
for the magnetic mass and z ~ 5 for the electric mass). In this region fits with 6 = and 
b 7^ yield mutually consistent results within statistical errors. From fits with b ^ we 
find best fits with b < 0, in accordance with the behaviour of local masses discussed above. 
Although these good quality fits with b < start at shorter distances, the magnitude of b 
is not well determined within our present statistical accuracy. In the following we thus will 
quote results from fits with 6 = 0. 

4.3 Numerical Results for the Propagators of the 3d SU(2) Higgs 
Model 

One of the aim of the numerical simulations is to clarify how far the screening mass can be 
determined in the framework of the effective 3d theory and to find the most suitable choice 
of h(x, 0). 

As it was already noticed before it is expected that the symmetric phase of the 3d model 
corresponds to the physics of the high temperature 4d SU(2) theory. Nevetherless let us 
investigate the propagators in the broken phase by fixing the parameter h to its value given by 
perturbative dimensional reduction h p Ad (x). In the broken phase the simulations were done for 
two sets of parameters: (3 = 16, x = 0.03, h = -0.2181 and (3 = 8 ,x = 0.09, h = -0.5159. 
The propagators obtained by us in the broken phase show a behaviour which is very different 
from that in the symmetric phase and that in the 4d case studied in Refs. [ESI HD]]. The 



magnetic mass extracted from the gauge field propagators is 0. 104(20)gf for the first set of 
parameters and 0.094(8)^1 for the second set of parameters. It thus is a factor 4 to 5 smaller 
than the corresponding 4d result. Moreover, the propagator of the A field does not seem 
to show a simple exponential behaviour, this fact actually is in qualitative agreement with 
the results of analytic calculations of Ref . H3[ . Taken together these facts suggest that the 



broken phase does not correspond to the high T physics of the 4d system. 

Now let us discuss the procedure of choosing the parameter h in the symmetric phase. 
The first step of this procedure is to chose for each value of a; a number of trial values for h 
for some fixed /3 and measuring the propagators for these parameters. The value of j3 was 
chosen to be (3 = 16. The second step might be an extrapolation of the screening masses 
measured for different values of h for each x to the value of the screening masses obtained 
in 4d simulation and determination of the corresponding value of h. The guiding principle 
in choosing the trial values of h was the fact that these values should be close to the value 
corresponding to the transition. This is based on observation that the perturbative line of 
4d physics lies close to the transition line. We suspect the same is true for the physically 
relevant h values with the exception that they are located in the symmetric phase. Therefore 



12 The definition of the local masses are given in Appendix B.2 where also futher details about the fitting 
procedure can be found 
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first we have to determine the transition line h c (x). The transition line as function of x in 
the infinite volume limit was found in |34| in terms of the renormalized mass parameter y 
(see appendix B.l for definition of y). The transition line in terms of y turns out to be 
independent of (3. Then using Eqs. (|B.39| ) and (|B.38| ) one can calculate h c (x). The usage 
of the infinite volume result for the transition line seems to be justified because most of our 
simulation were done on a 32 2 x 64 lattice. The two sets of h(x) values, which appear on 
Figure 10, were chosen so that the renormalized mass parameter y always stays 10% and 
25% away from the transition line. The values of the parameters in the symmetric phase 
are shown in Table 1 where the two sets of h values are denoted as (I) and (II) and also the 
values of h c corresponding to the transition line are given. Theses trial values together with 
perturbative line of 4d physics are also depicted in Figure 10. To check the suitability of trial 
values shown in Table 1 simulation for x = 0.07 and value of h = —0.2179 which lies deeply 
in the symmetric phase was done. The corresponding value of the Debye mass was found 
to be 2.41(11). This value should be compared with the value of the Debye mass obtained 
mo/T = 1.85 obtained in 4d case jEJ for T/T c = 12.57 which corresponds to x — 0.07. 
Clearly this choice of h is not suitable and one has to choose smaller h values to reproduce 
4d data. The analysis described above motivated our choice of trial values for h at other 
values of x as given in Table 1. 



Temperature scale 


h 


X 


T/T c 


I 


II 


transition 


0.09 


4.433 


- 0.2652 


- 0.2622 


- 0.2672(4) 


0.07 


12.57 


- 0.2528 


- 0.2490 


- 0.2553(5) 


0.05 


86.36 


- 0.2365 


- 0.2314 


- 0.2399(6) 


0.03 


8761 


- 0.2085 


- 0.2006 


-0.2138(9) 



Table 1: The two sets of the bare mass squared used in the simulation and those which 

correspond to the transition line for {3 = 16 

Yet another possiblity to be discussed is the simulation in the metastable symmetric 
phase along the line of perturbative 4d physics. As was already pointed out due to the fact 
that the transition line is strongly first order for all physically relevant values of x metastable 
symmetric branch exists below the transition line (see Figure 10). In fact the perturbative 
line of 4d physics lies in the region where this metastable branch exists. If one prepares the 
system to be in this metastable phase by starting a Monte-Carlo simulation with disordered 
random configuration it will not tunel to the other broken phase corresponding to the broken 
phase 0. Simulations in the metastable phase were done for parameters shown in Table 1. 
The temperature dependence of the screening masses obtained in the symmetric phase for 
two sets of the parameters, as well for the perturbative line of 4d physics h p Ad (x). in the 
metastable phase is shown in Figure 11. Also given there is the result of the 4d simulations 
49], m 2 D /T 2 = Ag 2 {T), with A = 1.70(2) for the electric mass and m T /T = Cg 2 (T), with 
C = 0.456(6) for the magnetic mass. As one can from Figure 11 the magnetic mass agrees 
quite well with the 4d data for both sets of h values given in Table 1 and it is not sensitive to 



13 The tuneling between two branches is supressed by a factor which is the exponential of the volume for 
all standard Monte-Carlo updating algorithms. 
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-0.28 



0.03 0.04 0.05 0.06 0.07 0.08 0.09 



Figure 10: The phase diagram of the 3d SU(2) adjoint Higgs model and bare mass parameters 
h used in our analysis (squares) for (3 = 16: perturbative line (a) and two sets / (b) and II 
(c). For a discussion of their choice see text. The thick solid line is the transition line. The 
dashed line (d) indicates the values of h below which no metastable branch exists and was 
estimated based on results of Ref. 



34 




10000 



T/T r 



Figure 11: The screening masses in units of the temperature. Shown are the Debye mass 
m D for the first (filled circles) and the second (open circles) set of h, and the magnetic mass 
rriT for the first (filled squares) and the second (open squares) set of h. Also shown there are 
the values of the Debye mass measured in metastable phase along the perturbative line of 
4d physics (open triangles). The line (a) and line (b) represent the fit for the temperature 
dependence of the Debye and the magnetic mass from 4d simulations from J|9|. Some data 
points at the temperature T 
scale for better visualization. 



90T C and 



9000T C have been shifted in the temperature 
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the actual value of h. The magnetic mass calculated for the perturbative line of 4d physics 
yields roughly the same result. The electric mass shows more pronounced dependence on 
h for temperatures smaller than 15T C and as one can see from Figure 11 only the values of 
the electric mass calculated for h p Ad are consistent with the 4d data. For larger temperatures 
the dependence of the electric mass on the actual value of h weakens but the second set of 
h values, which lies more deeply in the symmetric phase, seems to overestimate the value of 
m D . 

Based on these observations we can conclude: the the interpolation procedure for finding 
the true value of h^d(x) is not necessary and one can assume that the correct choice is 
h? Ad (x). This choice for the line of 4d physics reproduces the screening masses with about 5% 
accuracy in the entire temperature range. This is exactly the precision which is reasonable to 
assume for the 3d effective theory. Nevertheless it is important to notice that an accepatable 
description of the screening masses can be also achived with the choice of h corresponding 
to the symmetric phase. 



4.4 Magnetic Mass in 3d SU(2) Pure Gauge Theory 



The magnetic mass found in the previous section seems to scale with the 3d gauge coupling 
g\ ~ g 2 (T)T. A similar behaviour was found for the square root of the spatial string tension 
in the deconfined phase of the finite temperature SU(2) gauge theory ||91|| . Moreover, the 
value of the spatial string tension is tJo~ s = 0.369(14)g 2 (T)T, which is very close the the value 
of the string tension calculated in the 3d pure SU(2) gauge theory: y/d^ = 0.3340(25)gf [j92|| . 
This behaviour can be understood for very high temperature (small coupling) because in this 
case, as it was discussed in subsection 2.3, the heavy A Q field with mass ~ gT can also be 
integrated out and the IR behaviour of high temperature SU(N) gauge theory is described 
by a 3d pure gauge theory in which the only dimensionfull scale is g\ ~ g 2 T. Unfortunately 
for the realistic temperature range (which is also studied here ) the above argument fail to 
hold because the coupling is large. A non-perturbative study is therefore needed to establish 
the relation between the magnetic mass found in the finite temperature SU(2) theory and 
the magnetic mass of the 3d pure gauge theory. 

We have measured the Landau gauge propagators for the 3d SU(2) gauge theory and 
from its large distance behaviour extracted the magnetic mass. The results for different 
values of (3 are listed in Table 2. 



p 


m T /gl 


X 2 /d.o.f 


12.00 


0.48 ±0.036 


0.580 


16.00 


0.42 ± 0.070 


0.497 


20.00 


0.44 ±0.068 


1.439 



Table 2: The results of the fit for the magnetic mass in 3d pure gauge theory. Simulations 

for (3 = 20 were done on 40 2 x 96 lattice 



Using the data from Table 2 one finds tut = 0.46(3) c/f. This value is in good agreement 
with the 3d adjoint Higgs model result. The magnetic mass thus is rather insensitive to the 
dynamics of the A field. This finding is in accordance with the gap equation study of the 
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adjoint Higgs model described in the previous section. Gauge boson propagators were also 



studied in the 3d SU(2) Higgs model in Ref. f38|. It was found that in the symmetric phase 



the magnetic mass extracted from the propagator is insensitive to the values of the scalar 
couplings and its value is 0.35(l)gf. The magnetic mass in the 3d pure gauge theory was 
also measured there and its value was found to a bit larger but close 0.35(1) (?f (see Figure 
3 in Ref. p3|). This value is considerably smaller than the value found by us. This fact is 



probably due to small lattices used in Ref. [53] and to the different method for extracting 



the screening masses. These issues are discussed in appendix B.2. 

The decoupling of static magnetic sector is also seen in study of gauge invariant cor- 



relators. In Refs. f34] , |94[ the mass of the lowest lying gluball was measured in the 3d 
SU(2) adjoint Higgs model and its value was found rather insensitive to the values of scalar 
couplings and close to the corresponding value of the 3d pure SU(2) gauge theory. 

4.5 Gauge Invariant Correlators and the Additive Constituent 
Gluon Model 

In this subsection I will establish a simple relation between gauge invariant and gauge de- 
pendent correlators. Gauge invariant correlators for the SU(2) adjoint Higgs model were 



measured in Refs. [0, 0. The correlators of the following gauge invariant operators were 
studied: Tr A 2 , h, tj = TrA i^ and G = F^Ffi Q 

Let us first consider the correlation function of the scalar operator TrA^ whose large 
distance behaviour gives the mass of the A —A bound state: m(A ). The mass of this bound 



state is also expected to determine the exponential falloff of the Polyakov loop correlator [^4 
Let us discuss the relation of m(A ) and the Debye mass defined through the propagator. 
If the coupling constant is small enough then m(A ) ~ 2m£, Q . In our case this perturbative 



relation is not satisfied. However, following the analysis of Ref. |TB[ the mass of the A — A 
bound state can be represented by the sum of two constituent adjoint scalar masses, which 
is mo- 
Next we discuss the gauge invariant correlator contaning gauge field operators. In Ref. 
T8fl it was suggested that the screening masses extracted from these correlators can be related 



to the screening masses extracted from propagators using the notion of the constituent gluon. 
According to this idea one can associate with each covariant derivative in the corresponding 
operator a constituent gluon with the mass equal to the propagator pole mass my. It was 
shown in Ref. |18|] that using this simple idea the spectrum of 3d SU(2) Higgs model can 
be qualitatively well interpreted. The mass of the lowest lying glueball state was measured 
in Ref. |94| and was found to be rather independent of the couplings h and x of the A 
field, the value of this mass is wiq ~ 1.67(2) In terms of the constituent model this state 



can be viewed as a loose bound state of four constituent gluons |18] with a constituent mass 
equal to my, thus vtiq ~ 4my. If one uses the value of magnetic mass in pure gauge theory 
found in Monte-Carlo simulations one finds mg = 1.84(12) t^f. For the value of the magnetic 



mass given by Nair [|8T| one gets mg ~ 1.28(73. 



14 In Rcf.Q] the correlation function of the oparator hi — e^Tr A^Fjk was studied instead of hij. But 
the masses extracted from this correlation function are close to those extarcted from 
15 Earlier measurements of Ref. [B4| overestimated the mass of this state. 
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An attempt for a gauge invariant definition of the constituent gluon mass was given in 
Ref. [ 31 1 . The approach is based on the expactation value of the following operator 



G Fijkl {x,y) =< F«(x)U ab (x,y)F b kl (y) >, (83) 

where 

U ab (x,y) = 2Tr T a U(x,y)T b rf(x,y), 

U(x, y) = P exp (ig J" dx^T^ (84) 

with T a being SU(N) generators. For large separation one has Gp ~ exp(— M s \x — y\). 
According to the additive constituent gluon model M s = Imp. The numerical value was 
found to be M s ~ 0.73gg. Based on the additive constituent gluon model and the numerical 
value of the magnetic mass in pure SU(2) gauge theory one would expect M s ~ 0.92(6)#f. 
On the other hand with the value of the magnetic mass obtained by Nair [|8T| one expects 
M s = 0Mgi 

Finally we consider the correlation function of hij. This correlation function was used 



for non-perturbative definition of the Debye mass |3(J [H], |32|. In terms of the additive 
constituent gluon model the mass extracted from this correlation function is m/, = mp> + 
2mp. The masses predicted by the constituent model compared with the results of direct 
measurements are shown in Table 3. 



parameters 


m(A )/gl 


m h /gi 


(3 x h 


measured constituent model 


measured constituent model 


9 0.10 - 0.4764 
16 0.09 - 0.2622 
16 0.05 - 0.2314 
9 0.04 - 0.2883 
32 0.04 - 0.1247 
24 0.03 - 0.1475 


1.01(2) a 1.75(18) 
1.54(15) 1.78(6) 
2.28(20) 2.38(16) 
2.41(2) a 2.74(2) 
2.20(ll) b 2.42(1) 
3.03(65) 3.28(30) 


1.39(6) a 1.79(15) [ 1.51(1) ] 

2.16(20) a 2.29(7) [ 2.01(1) ] 
2.42(15) 6 2.12(7) [ 1.85(1) ] 



Table 3: The masses of the adjoint scalar and heavy-light bound states in units of g\ 
compared with the predictions of the aditive constituent gluon model. Some bound state 
masses were taken from Ref. |M| (a) and Ref. ( D ), other values were taken from Ref. 
BTfl . In square brackets we indicate the predictions of the constituent gluon model based 
the value of the magnetic mass estimated by Nair 



As one can see from the Table except for the largest value of x which corresponds to the 
temperature T ~ 3T C the simple additive model can describe the spectrum of gauge invariant 
screening masses with the accuracy of 15 — 20%. This is the accuracy which is reasonable 
to expect from such simple bound state picture which does not take into account binding 
energy. 

The success of the additive model in describing the mass of the A — A bound state 
has important physical consequence concerning the physics of the screening defined from the 
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Polyakov loop correlator. In section 2.4 it was argued that the true large distance behaviour 
of the Polyakov loop correlator is determined by the lightest 3d gluball state. However, as 
one can see from Table 3 for x < 0.5 the the masses extracted from the Tr Aq correlator 
are larger than the 3d gluball mass and can be well described by the additive model. This 
suggests that the dominant contribution to the Polyakov loop correlator comes from the 
exchange of two constituent A field in a similar manner as it happens in the leading order 
of perturbation theory. 



4.6 Gauge Dependence of the Screening Masses 



Before closing this section we briefly want to address the question of a possible gauge de- 
pendence of the screening masses extracted by us. The pole of the propagator was proven 
to be gauge invariant in perturbation theory (E2l |9~5[]. The statment on the propagator pole 



masses extracted from lattice simulation is less clear. Here one faces the numerical problem 
to isolate the asymptotic large distance behaviour of the correlation function from possible 
short distance (powerlike) corrections which are gauge dependent p8| . In Ref. |89| quark- 
propagators were studied in axial, Coulomb and Landau gauges and the effective masses 
extracted at quite short physical distances were found to depend on the gauge. In Ref. |38) 
the quark and gluon propagators were investigated in so-called A-gauges. Masses extracted 
from the propagators using exponential fits like Eq. fl32|) also show a mild dependence on 
the gauge parameter. To study the gauge dependence of our results, following, 
introduced A-gauges which in our case are defined by the condition 



we have 



Xd 3 A 3 + d 2 A 2 + dxAi = 0. 



(85) 



The propagators have been measured in A gauges and the masses were extracted from the 
propagators using the functional form given in Eq. fl32|). 

Some preliminary results on the gauge dependence of the screening masses were already 
published in 



However the conclusion presented there are somewhat different from the 
present situation due to the low statistics of the first analysis. 

For the numerical analysis the following values of the parameter A have been chosen : 
A = 0.5, 1.0, 2.0, 8.0. Simulations were performed on a 32 2 x 96 lattice at x = 0.03 and for 
two values of (3 and h. Values of the Debye mass measured for different A and two values of 
/3 are summarized in Table 4. 



A 


0.5 


1.0 


2.0 


8.0 


@ = 16, h = -0.2085 
= 24, h = -0.1510 


1.28(6) 
1.26(7) 


1.21(9) 
1.31(9) 


1.15(5) 
1.20(5) 


1.08(6) 
1.11(3) 



Table 4: The Debye mass in units of the temperature for x 

different A parameters. 



0.03 and calculated for 



As one can see from the results the Debye mass shows a weak dependence on A. Its value 
is decreasing as A increases but for A > 2 it seems to approach a constant value which is 
independent of A. This dependence on A is similar to those which was found in Ref. [38]. 
In Figure |T2| we show the local electric masses for (3 = 24 and h = —0.1510 for the largest 
(A = 8) and the smallest (A = 0.5) values of the A parameter. 
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Figure 12: The local electric masses for f3 — 16, x — 0.03 and h = —0.1510 in units of g\ 
measured on 32 2 x 96 lattice for different values of the gauge parameter A. Shown are the 
local electric mass for A = 0.5 (a) and for A = 8.0 (b). 



N 




Figure 13: The local magnetic masses for /3 — 16, x — 0.03 and y = 0.9279 in units of g\ 
measured on 32 2 x 96 lattice for different values of the gauge parameter A. Shown are the 
local magnetic mass for A = 0.5 (a), for A = 1.0 (b) and for A = 2.0 (c). 
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The local magnetic masses measured at (3 — 16, x — 0.03 and h = —0.2085 for A = 
0.5, 1.0, 2.0 are shown in Figure 13. Similary to the electric mass the magnetic mass also 
show a mild dependence on A. 

Futher studies are required to clarify the source of the non-perturbative gauge dependence 
of the screening masses. 
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5 Coupled Gap Equation for SU(2) Higgs Model 



As it was already mentioned in the Introduction, considerable progress has been achived in 
understanding the thermodynamics of the electroweak phase transition in the past 5 years. 
For the electroweak theory the separation of different mass scales holds (2irT 3> gT ^> g 2 T)f^\ 
and the superheavy modes (i.e. the non-zero Matsubara modes with typical mass ~ 2ttT) 
together with the heavy A field (with a mass ~ gT) are integrated. The thermodynamics is 



described by an effective theory, the 3d SU(2) Higgs model |T0|, [TT|, [12|, [40[| . The properties 



of the phase transition and the screening masses were studied in great detail using lattice 
Monte-Carlo simulations of the reduced model |19|, |93|. The application of the gap equation 
approach [ |68| , [T0|] led to similar results. Both approaches predict that the line of first order 
transitions ends for some Higgs mass m# = m° H and the value of the critical mass m c H was 
found to be the same within 20%. 

Recent 4d Monte-Carlo simulations of the finite temperature SU(2) Higgs model [97, |98| 
provide good non-perturbative tests for the validity of dimensional reduction. A detailed 



discussion of relating 4d and 3d results was published very recently in Ref. |99[ . 

The purpose of the investigations presented in this section is twofold. The first aim is to 
provide some non-perturbative evidence for the decoupling of the A field from the gauge + 
Higgs dynamics in the vicinity of the phase transition. The coupled set of gap equations for 
the screening masses will be derived and solved numerically in the 3d fundamental + adjoint 
Higgs model. This model emerges when the non-static modes are integrated out in the full 
finite temperature Higgs system. Its predictions for the screening masses will be compared 
with those obtained by Buchmiiller and Philipsen (BP) [[F(J in the 3d Higgs model using the 
same technique. This model is obtained from the 3d fundamental + adjoint Higgs model by 
integrating out the heavy Aq field. The main result of this comparison is a proposition for 
a non-perturbative and non-linear mapping between the two models ensuring quantitative 
agreement between the screening masses in a wide temperature range on both sides of the 
transition. This high quality evidence for the decoupling of the A field at the actual finite 
mass ratios presumes, however, the knowledge of the "exact" value of the Debye screening 
mass, since for the proposed mapping its non-perturbatively determined value turns out to 
be essential. 

The second aim of the present investigations is to examine the symmetric phase in more 
detail. There the Higgs and the Debye screening masses are both of the same order of mag- 
nitude ~ gT and thus in that regime there is no a priori reason for the A field to decouple. 
This circumstance makes the quantitative relation of the screening masses calculated in the 
3d fundamental + adjoint Higgs model particularly interesting in the high-T phase. Here we 
are going to apply two different resummation techniques and check to what extent persists a 
non-perturbative mass hierarchy in this part of the spectra. Investigations presented in this 
Section are based on Ref. ||100| . 



16 The renormalized gauge coupling for the electroweak theory has a moderate value g ~ 0.66. In the 
following discussion we will ignore the presence of the U(l) sector and fermions because they are not essential 
from the point of view of basic properties of the electroweak phase transition. 
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5.1 The Extended Gap Equations 



The Lagrangian of the three dimensional SU(2) fundamental + adjoint Higgs model is J70|, |TT 



L 3D = Tr[-F ij F ij + (D i ^) + (D i ^)+f/^ + ^ + 2X(^ + ^) 2 ) 

+\{D l A f + \AA + ^{Alf + 2ci 2 Tr$+$, (86) 



where 



$ = 1(^1 + ^), A$ = (^-W)$, W i =^7tf i . (87) 

The relations between the parameters of the 3d theory and those of the 4d theory are derived 
at 1-loop level perturbatively fTTfl: 



c = yid T > rfb = yidT 2 ' V 2 = {j&sld + |^4d) T 2 - \[i\ d . (88) 

If the integration over the A adjoint Higgs field is performed we obtain the model investi- 
gated in Jf(J with parameters g, A, p,. These couplings of the reduced theory are related to 
the parameters of the 3d fundamental + adjoint Higgs theory through the following relations: 

-2 2 fi 9 2 \ T \ 3C 2 _ 2 2 3cfI D 

g=9 \ l -2^- D )^ A = A "W "^""aT (89) 

In particular, we note that the /2 scale serves as the temperature scale of the fully reduced 
system, while \x is the scale for the system containing both the fundamental and the adjoint 
scalars. The two are related perturbatively by a constant shift. 

In order to perform the actual calculations in the broken phase it is necessary to shift 
the Higgs field, a — > v + a'. After this shift and the gauge-fixing (the gauge fixing parameter 
is denoted by £) the Lagrangian including the ghost terms assumes the form 



2 

+ 9 -va'Wl + |w> M • (7r<V " °'d^) + |(H> M x tt) ■ dfi 



+ jW 2 (a' 2 + Tf 2 ) + \va\a' 2 + tt 2 ) + -^a' 2 + vr 2 



+-(AA) +-m 2 m A +— (Aq) +2cva'A + cA Q (a' 2 + vf 2 



(90) 
(91) 
(92) 
(93) 
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q 2 q 2 11 

+gd^c* ■ (Wp x c) + i^-va'c*c + i^-vc* ■ (vf x c) + ~/iV + -Xv 4 



(94) 
(95) 

(96) 

+ i(/i 2 + A^ 2 )(a' 2 + vf 2 ) + v(fi 2 + \v 2 )a' , (97) 

where the following notations were introduced for the tree-level masses: m^ = \g 2 v 2 (the 
vector boson mass), Mq = /i 2 + 3At> 2 (the Higgs mass) and m 2 m — p? D + 2ct> 2 (the Debye 



mass). The last two terms of fl97|) arise from the Higgs potential after the shift in the Higgs 



field a. For /i 2 < 0, they vanish if one expands around the classical minimum v 2 = — /i 2 /A. 



These terms have to be kept [70], when the equation for the vacuum expectation value is to 
be deduced. 

In order to obtain the coupled gap equations one replaces the tree-level masses by the 
exact masses 

m 2 TQ ^m 2 T + 5m 2 T , M 2 —> M 2 + 5M 2 , m 2 DQ -> m 2 D + 5m 2 D , (98) 

and treats the differences dm? = m^ — m 2 ^, 5M 2 = Mq — M 2 , 5m 2 D = m 2 m — m 2 D as 
counterterms. The exact Goldstone and ghost masses are both equal to ^tut, where is 
the exact gauge boson mass. The gauge invariance of the self-energies of the Higgs and gauge 
bosons is ensured by introducing appropriate vertex resummations. Their explicit formulae 
can be found in [ffO| . In the present extended model, a resummation of the Higgs-A vertex 



would be also necessary if the gauge invariance of the Aq self-energy is to be ensured. Then 
the only source of the gauge dependence which would remain is the equation for the vacuum 
expectation value v. 

All these resummations are equivalent to work with the following gauge invariant La- 
grangian: 

Lf = \P 13 F 10 + Tr ((A*) + A$ - \M 2 ^ +\{m 2 D - ^)A 2 



o 2 M 2 

+ ^-^-Tr($ + $) 2 + 2cA 2 Tr$ + $. (99) 
In this Lagrangian one shifts the Higgs field around its classical minimum a — » a' + and 



adds the corresponding gauge fixing and ghost terms |70| . Shortly, we shall argue that the 
A -Higgs vertex resummation arising from the replacement of v by 2rriT/g when the scalar 
field is shifted in the last term of the above Lagrangian destroys the mass-hiearchy between 
the heavy A and the light gauge and Higgs fields. Therefore we have to give up the full 
gauge independence of the resummation scheme. The numerical solution to be presented 
below shows that the gauge dependence of the Aq — $ vertex in our resummation scheme 
introduces only a minor additional gauge dependence beyond that of the equation for the 
vacuum expectation value |70f appearing below in Eq. ( |106|) . 
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The coupled set of gap equations is constructed from that of Ref. JF(J by adding the 
contributions due to the presence of the adjoint Higgs field. The self-energy contributions 
for the 3d adjoint Higgs model were already calculated in section 3. Below we list only the 
additional contributions to the self-energies, which all contain at least one Aq — $ vertex (the 
corresponding diagrams are listed in the Appendix B). Let us emphasize once again that no 
resummation of the A — $ vertex was applied. 

The additional contribution to the self-energy of the A field coming from Higgs, Gold- 
stone, gauge and ghost fields (diagrams a-i) is 



xy\H,G,gh 
0LL A 



Acv 2 (/i 2 + Xv 2 ) 3cgv 



M 2 



71 



m 



r 



M 

4^ + M 2 



cM 3v^ 



Ac 2 v 2 



71 



3 mi) 1 p 

_ ^tt77 arctan — 

2 M 2 p m D + M 



(100) 



There is also an additional contribution to the gauge boson self-energy coming from the 
adjoint Higgs field (diagram m): 



511% (p, mr, M, mo) 



3cg tutv 



-m D . 



2n M 2 

The contribution of A to the Higgs self-energy (diagrams j-1) is the following: 

3m dc 6c 2 v 2 1 



sag fa 



m T ,m D 



2tt 



-- arctan • 



71 p 



P_ 

2m 



D 



9gcv 

+ -: rn D . 

Anmr 



(101) 



(102) 



Making use also of the pieces of the self-energies calculated in section 3 and also in [70] we 
write down a set of coupled on-shell gap equations for the screening masses of the magnetic 
gauge bosons, fundamental Higgs and adjoint A fields in the form: 



m q 



M 



m 



D 



U T (p 2 = -m 2 T , m, M) + <ffl^° (p 2 = -ml, m D ) 

(p 2 = -m 2 T , m T , M, m D ) (103) 

S(p 2 = -M 2 ,m T ,M) + 5Tlp(p 2 = -M 2 ,m T ,m D ), (104) 

n o(p 2 = m Di m T , m D ) + 5n%f' gh (p 2 = - m Di m T, M, m D ), (105) 



where and £ are defined by Eqs. (17), (18) of Ref. fF(J. 5U^° and n 00 are defined by 
Eqs. (§7[) and (gj). 

If on the right hand side of the third equation one inserts the tree level masses, the 
next-to-leading order result of Ref. | |101| | is recovered for the Debye mass in the SU(2) Higgs 
model. 

The equation for the vacuum expectation value completes the set of the above three 
equations: 



v (/i 2 + Xv 2 ) 



16tt 



g (^m 2 T + v^M 2 + 



M 3 \ 3c 
+ —vm D . 

mr J 27T 



It is important to notice that this equation can be rewritten as 

3 



^9 Um% + VIM 2 + 
lo7r \ mr ) 



(106) 



(107) 
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Figure 14: The Higgs mass in units of g 2 as function of /J?/g calculated at X/g 2 = 1/8 using 
the gauge-invariant Aq — <1> vertex resummation version of the gap equations. Shown are the Higgs 
mass derived in the full static theory in the Landau- (a) and in the Feynman-gauge (b) and the 
Higgs mass, and in the ^-reduced theory in the Landau-gauge (c) and in the Feynman-gauge (d). 
The // 2 -shift indicated by eq,(|89|) was applied. 



2 _ 2 



with 

tits = ^r mD - ( 108 ) 

This equation is formally identical to the equation of BP for the vacuum expectation value 
70fl . On the basis of this observation, we expect that the main effect of the Aq integration 



is the above shift in the /z 2 -scale. Since mj) is itself a non-trivial function of \x this non- 
perturbative mapping is also nonlinear. 

A very similar set of equations could be derived for the case of the gauge invariant 
resummation of the Aq — $ vertex. For instance, one would write in the last term on the 
right hand side of ( |106| ) 2rriT/g on the place of v , which would prevent the absorption of this 
term into a redefinition of the temperature scale. This was the reason when we solved the 
gap equations first for the case of gauge invariant Aq — $ vertex resummation, that we have 
found substantial corrections to the BP results, beyond the shift of the /x 2 -scale indicated in 
eq. (|89f) . This situation is illustrated in Figure |14[ 

This makes it clear that the decoupling of the y4 -field for the physical value of the pa- 
rameter c would not work, the hierarchy of the screening masses would be destroyed under 
such vertex resummation scheme. Since our goal is to propose such a mapping onto the fun- 
damental Higgs-model, which already for finite y4 -mass shows evidence for the asymptotic 
Appelquist-Carazzone theorem we have renounced from the gauge invariant resumma- 
tion of the Aq — $ vertex. 



5.2 Numerical Results 

The main goal of the present investigation is to analyze non-perturbatively in quantitative 
terms the decoupling of the Aq field. We are interested in those effects which go beyond the 
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Figure 15: The Higgs boson masses at X/g 2 = 1/8 (crossover region) in units of g 2 as function 
of /J?/g 4 in the 3d fundamental + adjoint Higgs model and in the 3d SU(2) Higgs (.Ao-reduced) 
theory. Shown are the Higgs mass in the full static theory in the £ = (Landau) gauge (a) and in 
the £ = 1 (Feynman) gauge (b), and the Higgs boson mass in the Ao-reduced theory in the £ = 
gauge (c) and in the the £ = 1 gauge (d). 



perturbative mapping implied by eq. (p9|) . Therefore, we will first compare the predictions 
for the Higgs and gauge boson masses from the coupled gap equations ( |103| )- (|106|) of the 
3d fundamental + adjoint Higgs model with those obtained in the ^-reduced theory, the 
3d Higgs model |f70| . The corresponding Higgs masses are shown in Figure [15] using two 
different gauges. The results obtained in the ^-reduced theory are displayed after the shift 
required by eq.fl89|) is performed. The difference between the full and the reduced theory is 
still visible in the vicinity of the crossover. In this region the relative difference between the 
prediction of the full and the reduced theory is about 20%. 

Our proposal to resolve this relatively large deviation is to introduce a more complicated 
relationship between the couplings. Having gained intuition from eq.( |108| ), we have plotted 
the mass-predictions for the Higgs-field derived from our full set of equations against the 
results of BP calculated for couplings taken from (|89|) with a replacement /i^ — > m^: 



The non-trivial nature of this replacement becomes clear from Figure 16 where the fi 2 - 
dependence of mjj is displayed. Clearly, its non-trivial /independence is most expressed 
in the neighbourhood of the phase transformation (crossover) point ^ 2 /g 4 G (0.1 — 0.2). 
The application of this mapping to the data obtained from the model containing both the 
fundamental and the adjoint representation leads to a perfect agreement of the two data sets 
for large values of X/g 2 . For smaller values of \/g 2 (1/32,1/64) the mapping ( |109|) works 
very well in the symmetric phase, but in the broken phase (|8§D seems to be the better choice. 



We suspect, that the tree level piece in mu arising from the Higgs-effect, should not be 
included into the correction of (^), since it is itself a tree-level effect. Therefore we propose 
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Figure 16: The fi 2 dependence of the Debye mass for ^ = 1/8. 



the following replacement in (|109|) : 



ni 



D 



2cv 2 . 



110) 



In Figure [H] it is obvious that a very good agreement could be obtained with this mapping 
between the Higgs mass predictions of the one-loop gap equations of the full static and the 
yl -reduced theory for X/g 2 = 1/32. The quality of the agreement on both sides of the phase 
transition is good, signalling that the influence of the "mini-Higgs" effect in the symmetric 
phase is negligible. Therefore it is not surprising that for X/g 2 = 1/8 the same quality of 
agreement is obtained like before. 

It is important to notice that there is a strong gauge parameter dependence in the 
symmetric phase and in the vicinity of the crossover. The variations due to the change in 
the gauge are equal in the full and in the reduced theory, which indicates that the additional 
gauge dependence, introduced by the gauge non-invariant resummation of the A field is 
negligible. The mapping ( |110| ) performs equally well in Landau- and Feynman-gauges. 

Other quantities which are worth of considering for the comparison of the full 3d and 
the reduced theories are X c /g 2 , the endpoint of the first order transition line and fi+/g 2 , the 
mass parameter above which the broken phase is no longer metastable. The values of fJ^/g 4 
for different scalar couplings and different gauges in the full and in the reduced theory are 
summarized in Table 1. Here the mapping (|109|) could be implemented only by extrapolating 
from smaller fi 2 /g 4 , since the end-points of metastability do not correspond to each other, 
and in some cases mo could not be determined from the gap equations. Also here for larger 
values of X/g 2 the application of ( |109| ) led to an improved agreement between the end-point 
l^+/g A values, while for X/g 2 = 1/48, 1/64 the mapping fl89|) works better. In the table we 
have displayed fi 2 ^/ ' g A values of the Ao-reduced theory shifted perturbatively and with help 
of the best performing non-perturbative mapping ( |110|) . For both gauges the latter agrees 
with the fi^/g 4 - values of the full static theory very well. 
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Figure 17: The Higgs boson mass in units of g 2 as function of fJ, 2 /g 4 calculated at X/g 2 = 1/32 
in the Landau gauge in the full static theory and in the Ao-reduced theory. Shown are the Higgs 
mass in the reduced theory obtained by perturbative reduction (a), in the reduced theory obtained 
by non-perurbative matching (cf. eqs. ( |109| ), ( |11C| )) (b) and in full static theory (c). 



A/s 2 


A 


B 


C 


£=0 £=1 


C=0 £=1 


£ = t = \ 


1/32 
1/48 
1/64 


0.1516 0.1423 
0.1647 0.1558 
0.1841 0.1750 


0.1426 0.1341 
0.1627 0.1541 
0.1881 0.1808 


0.1499 0.1405 
0.1637 0.1546 
0.1875 0.1792 



Table. 1: Values of fi+/g in the full static theory (A), in the perturbatively reduced theory (B) 
and in the reduced theory obtained using non-perturbative matching described in the text (C). 
Calculations were done in the Landau (£ = 0) and in the Feynman (£ = 1) gauges. 

The endpoint of the 1 st order line in the Landau gauge in the 3d Higgs theory was found 
at X c /g 2 = 0.058. The corresponding critical scalar coupling in the full 3d theory is within 
the 1% range. In Feynman gauge we find X c /g 2 = 0.078 for the Ao-reduced theory and 
the corresponding value for the full 3d theory lies again very close to it. Thus the A field 
has almost no effect on the position of the endpoint. The strong gauge dependence of A c 
indicates, however, that higher order corrections to this quantity are important. 

The gauge dependence of the screening masses is even more pronounced deep in the 
symmetric phase (fi 2 /g A > 0.3). For example the value of the gauge boson mass is roughly 
0.28g 2 in the symmetric phase for the Landau gauge. The corresponding value in the Feyn- 
man gauge is about 0.22g 2 . The gauge dependence of the gauge boson mass is somewhat 



weaker at the 2-loop level ||80| . It should be also noticed that the gauge boson mass depends 
weakly on the parameters of the scalar sector (/i, /x D , X, Xa)- This fact was also noticed in 
previous investigations |70], [76 . 
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5.3 Screening Masses in the Symmetric Phase with a Gauge In- 
variant Resummation Scheme 

The main motivation for the present investigation was to gain insight into the decoupling of 
the dynamics of the fundamental and the adjoint Higgs fields. The degree of the decoupling 
is expected to depend on the mass ratio of the fundamental and adjoint Higgs fields. In the 
symmetric phase both masses are of the same order in magnitude (eg. ~ gT) . Therefore the 
hierachy of the A and Higgs masses can only be present due to numerical prefactors. The 
persistence of the perturbatively calculated ratio should be checked in any non-perturbative 
approach. 

As we have seen in the previous section the gauge dependence in the symmetric phase 
is too strong in the applied schemes to give a stable estimate for the mass ratio of the 
fundamental and the adjoint Higgs fields. A reliable non-perturbative estimate for the Higgs 
mass deep in the symmetric phase (defined through the pole of the propagator) is even more 
interesting because it was not measured so far on lattice. Therefore, in this section we will 
investigate a coupled set of gap equations in the symmetric phase which is based on the 
gauge invariant resummation scheme of Alexanian and Nair (AN) W2\. In this approach 
one can avoid any vacuum expectation value for the Higgs field in the symmetric phase and 
because of this fact this approach is gauge invariant. 

To derive the gap equations in this scheme one has to add and substract the mass 
generating term defined by Eq. ( 53 ) and also add the corresponding gauge fixing term 
defined by Eq. (69). Among diagrams shown in Appendix B those which contain gauge 
boson propagators should be reevaluated. Straightforward calculations lead to the following 
equations: 



M 



m 



D 



2 

Cg 2 m T + g J^L(2f{m D /m T ) + f{M/m T ) 
+ T- l-,g 2 MF{M/m T ) - 6AM - 6cm D ) , 

4"7T \4 / 

fi 2 D + —(^2g 2 mDF(mD/iTLT) — 5Aa^d — 8cM^. 



Ill) 



where C = ^(21/4 In 3 — 1) |72[ and the following functions were introduced 

111 

f(z) = --z + (z 2 - -)arctanh— , 



F(z) = -1 - 



1 



(4z 



ln(l + 2z) 



112) 

;ii3) 



Let us first discuss the ratio of the A and the fundamental Higgs masses. In Figure [ISj 
this ratio is shown as calculated from the eqs. ( |1 1 1| ) and compared with the corresponding 
perturbative value. The /i interval in this plot corresponds to the temperature range relavant 
for the electroweak theory T < ITeV. We have also analyzed the //-dependence of the 
fundamental Higgs mass alone in the full static and in the y4 -reduced model. For [i 2 /g 4 in 
the interval (0.2 — 0.3) the result of the gauge invariant approach agrees fairly well with the 
masses obtained in the BP-scheme. In Figure [TP] the difference between the Higgs masses 
calculated from the coupled set of gap equations ( |1 1 1|) and the leading order perturbative 
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Figure 18: The ratio of the Debye and the fundamental Higgs masses for p- = 1/8 calculated from 
gap eqs. ( |H1| ) (a) and the leading order result (b). 
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Figure 19: The non-perturbative correction to the Higgs mass as the function of [i calculated from 
the full static (a) and from the Ag-reduced theory (b). 
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result (M ) is shown. As one can see, the non-perturbative correction to the Higgs mass is 
the largest for small ji and is decreasing as /i increases reaching the percent level for large 
enough ji. 

The relative difference between the full static and the A -reduced theory, however, is 
slowly increasing as \i increases and the hierarchy between the A and the Higgs masses 
becomes less pronounced as \x gets larger. The relative difference between the Higgs masses 
calculated in the full static and in the A -reduced theory varies between 20% for ji 2 j ' g A = 0.2 
and 35% for /i 2 /# 4 = 0.6. 

It is also important to notice that the A field is not sensitive to the dynamics of the 
Higgs field. In particular it turns out that m D depends weakly on ji and A in the symmetric 
phase and its value is close to the corresponding value calculated in 3d adjoint Higgs model. 
Let us notice that the magnetic mass in this resummation scheme also seems to be insensitive 
to the dynamics of scalars, therefore the magnetic and electric screening masses are close to 
their values determined in the pure SU(2) gauge model. 
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6 Conclusions 



In my Thesis I have studied the problem of screening in hot non-Abelian plasma. A de- 
tailed review of the phenomenon of chromoelectric screening was given. It was shown that 
the chromoelectric screening is a well understood phenonenon in the leading order of per- 
turbation theory and it is quite similar to the well known Debye screening phenomenon in 
non-relativistic electron plasma. Beyond the leading order, however, it turns out to be very 
complicated and at present it is even not fully clear how to define it properly. In my numeri- 
cal investigations Debye screening phenomenon was defined through the long distance (small 
momentum) behaviour of the static longitudinal propagator. It has been shown that such a 
definition is free of infrared singularities in higher orders of perturbation theory only if one 
introduces the so-called magnetic screening mass, a concept which was introduced long ago 
by Linde to cure infrared divergencies of finite temperature gauge theories [[T7|] . 

In recent years the problem of magnetic mass was actively studied using self-consistent 
resummation of perturbative series. The introduction of the magnetic mass allows a self- 
consistent determination of the Debye screening mass (the inverse screening length) through 
coupled gap equations, which was discussed in section 3. It was argued that such determi- 
nation so far is only possible in the framework of the dimensionally reduced effective theory. 
Thus one has to address the question how far the dimensionally reduced effective theory 
faithfully represents the high temperature limit of gauge theories. This question is far from 
being trivial since the separation of different mass scales 2nT ^> gT ^> g 2 T which is nec- 
essary for dimensional reduction is not satisfied for any realistic temperature range in the 
theory which corresponds to quantum chromodynamics because the gauge coupling constant 
g is large. This observation suggested that the electric (Debye) and the magnetic screening 
masses should be determined simultaneously from the coupled set of gap equation. In section 
3 where the screening masses were studied through coupled gap equation the precision of 
dimensional reduction was tested also by examining the effect of non-local operators and it 
turned out that their contribution is small thus justifying the approach based on dimensional 
reduction. 

In section 4 the screening masses were studied using lattice Monte-Carlo technique in the 
3d effective theory corresponding to the high temperature phase of SU(2) gauge theory, the 
3d adjoint Higgs model. Comparison of the screening masses obatined in these simulations 
with the corresponding masses obtained recently in the full 4d finite temperature theory 
49fl allows to test non-perturbatively and quantitatively the precision of the dimensionally 



reduced theory. It was shown that the predictions of the 3d effective theory are quite 
reasonable even at temperature ~ 4T C and become more and more precise as the temperature 
is increased. This fact gives prospects for studying the screening phenomena in the effective 
theory approach. The predictions of the gap equation approach were compared with the 
corresponding results. It turns out that the results obtained from coupled gap equations 
describe correctly the temperature dependence of the screening masses, however, so far there 
is no quantitaive agreement between the predictions of the coupled gap equations and the 
Monte-Carlo data. The Debye mass measured in lattice Monte-Carlo simulations is about 
25% larger than what is predicted from the coupled set of gap equations. 

In section 5 the screening masses of the finite temperature SU(2) Higgs model were stud- 
ied. Though the application of the dimensional reduction approach here seems to be more 
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justified (the corresponding gauge coupling constant is g ~ 0.66), the non-perturbative quan- 
titative test of the precision of the dimensional reduction is also important here. The study 
of the screening masses in this model through coupled gap equations tests the decoupling 
of the heavy Aq field (temporal component of the static gauge field). Detailed numerical 
analysis revealed high precision decoupling of the Aq field if in the mapping of the full static 
theory onto the v4 -reduced system its mass determined non-perturbatively is used. This 
observation also implies that for the realistic value of the mass perturbative integration 
over the heavy A Q field has limited precision due to the coupling of the heavy A Q field to the 
light magnetostatic sector. This coupling to the magnetostatic sector is taken into account in 
the proposed non-perturbative mapping procedure. It is very interesting that this coupling 
between the heavy A field and the light magnetostatic sector can described quite precisely 
by the magnetic mass obtained from the gap equations. 

From the above discussion it is clear that the concept of the magnetic screening mass 
is extremely usefull since it allows higher order definition of the electric screening mass, 
which quantitatively agrees with the corresponding lattice predictions and also the high 
quality mapping onto the effective theory mentioned above. The values of magnetic mass 
determined in different resummation schemes are quite close to each other. A recent 2-loop 
gap equation of the magnetic mass has shown that the value of the magnetic mass does not 
change dramatically relative to the 1-loop results and the corresponding value was found to 
be 0.34gf |8(| (which should be compared with the 1-loop predictions 0.28gf and 0.38gf). 



Moreover the non-perturbative analysis of the mass gap in 3d pure gauge theory yields a 
very similar value 0.32gf ||81| . All these facts give us confidence that the magnetic mass is a 
viable concept. However, in Refs. [75, 72] some inconsistencies in the 1-loop determination 
of the magnetic mass were found. In the numerical simulations the value of the magnetic 
mass was found to be ~ 0.46gf which is substantially larger than the predictions from the 
analytical approaches. Thus I think, further studies are required to clarify the phenomenon 
of non-Abelian magnetic mass generation. 
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Appendices 



A Calculation of Different Diagrams in SU(N) Adjoint 
Higgs Model 

The propagators of different fields can be read off the quadratic part of the Lagrangian. 
These exact propagators are listed below. The gauge boson propagator: 



D ^ k) ~ \ 5ij k 2 ) k 2 + m 2 T + k 2 k 2 + mY (A ' 1} 
where m^ = y^mr, the propagator for the adjoint scalar field A : 

k 2 + mjj 



and finally, the ghost propagator: 



m = jot (a.3) 



The Feynman diagrams contributing to the 2-point functions of the relevant fields are 
shown below. In these diagrams every line corresponds to a resummed propagator. The wavy 
line corresponds to the gauge bosons, the solid one to A and the dashed one to the ghost. 
The corresponding analytic expressions can be written in terms of the following standard 
integrals: 



I(k,r, s) 



I m (k,r,s) = 



d d p 1 
(27r) d (p + k) 2r p 2s 

k d-2(r +S ) r(r + s-f) r(f-s)r(f-r) 

(47r) d / 2 r(r)r(s) T(d-s-r) 
d d p p m 



(A.4) 



(27r) d (p + k) 2r p 



2s 



I mn (k,r, s) 



k d ~ 2 ^ T(r + s-l) r(f + I- g )r(f-r) 
( 4vr )d/2 T(r)T(s) T(d+l-s-r) ' 1 ' 

d d p p m p n 



(2ir) d (p + k) 2r p 2s 
k d-2(r+s) r p r (r + s - 1 - f ) r (f + f - s ) r (f + I - r) 
(47r) d / 2 12 T(r)T(s) T(d+2-s-r) 

r(r + 8 -i) r(f + 2- s )r(f-r) 



s 



mil 



r(r)r(s) T{d + 2 -s-r) 



k m k r 



(A.6) 



m s f d3 P 1 A 

J(k,m 1 ,m 2 ) = J ; ; ^ — 9 ; = — , (A.7) 



J m (k,m 1 ,m 2 ) 



if{{p + k) 2 + m 2 ){p 2 + m 2 ) 8tt' 

d 3 p Pm 

{2nf {{p + k) 2 + m 2 ) {p 2 + m 2 ) 
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J mn (k,m 1 ,m 2 ) = 



8^ 
d 3 p 



m l — 777 2 
k 2 



k 2 + to? — ml . 

-A 



2k 2 



(27r) 3 ((p + k) 2 + to 2 ) (p 2 + ml 



5' 



8tt 
8tt 



to x 4/c 2 to 2 + (/c 2 + to 2 — to 2 .) 2 
~2~ + 8k 2 A 



mi 
2k 2 



^k 2 m\ + 3 (k 2 + m\ — ml) 2 



8k 4 



1 (toi — to 2 ) (A; 2 + to 2 — ml) 



8n 



4A; 2 



k m k r ' 



where A = T arctan 

k 



m 1 +m,2 



j(m) 
l mn {m) 



d 3 p 



(2ir) 3 p 2 + to 2 

d 3 p PmPn 

(2n) 3 p 2 + m 2 



TO 
471"' 



m 2 j(m). 



These integrals were evaluated using dimensional regularization. 

The diagrams contributing to the 2-point functions of A^ are the following: 



a 



a 






(A.8) 
(A.9) 



(A.10) 



(A.H) 
(A.12) 



d e 

with the analytical contribution: 



n«7(*) 



= # 2 iV(5 a6 J mri (A:, TO G , 777.(5 ) + k m J n (k, m G , TO G )J , 



mn 



<7 2 A^™ -i(TO T ) + -O'K) 



fTO 2 + k 2 ) 2 



mn 



(J(k, to t , 0) + J(k, 0, to t )) 



2fc 2 (4to 2 + k 2 ) T/ , . / k 2 , 

+ % J -J(k,m T ,m T ) -2 (l + — ) j(m T ) 



TOq 



TOq 



(A.13) 
(A.14) 
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(mt + k 2 ) 2 T/J / k A k 2 \ 

4 J(k : m T ,0)+[—r- — -6) 



4m 4 



J(fc, m^, m^) 



+ 



+ 



^A+l^ + l) J (k, 0, %) + ^(mr) + ^I(k, 1, 1) 
A; 4 / k A k 2 



T 



2 \ 2 



lmn(jn T ) - [ 1 + — J I (Jmn(k, 0, m T ) + J mn (fc, m T , 0)) 



+ 



^rj(m T ) + -^kmlnik, 1, 1) - + + 3 ) k mJn{k, m T , 0) 



+ 1 



A' 1 



m 



f k A k 2 \ l-i 
k m J n (k, 0, m T ) + I — + 4— + 8 J k m J n (k, m T , m T ) \ 



-g 2 N5 ab i [5 n 



(m 2 + k 2 ) 2 



m, 



(J(k, 0, m T ) - m L , mr)) 



mf- 



k k 



Am\m\ 



m\m\ 



6k 2 m1 + 7m 4 -, + 1m\m\ — (m 2 L + A; 2 ) 2 j J(k, mj, mr) 

+m 2 L j{mT) — lm\.j[mj)j + [k A — 6k 2 m1 — 7m 4 -) J(k, 0, m r ) 
-A; 4 /(A;, 1, 1) + (mi + k 2 ) 2 J(k, m L , 0) 

m 2 T k m k n j{m L ) + A; 4 (J mn (k, m L , 0) - / mn (fc, 1, 1)) 



+ 



+ (k 2 + m|) (J mn (k, 0, m T ) - J mn (A;, m L , m T )) + m 2 T l mn {m L ) 
L ^-k m J n (k,m L ,m T ) ^-^k m I n (k, 1, 1) 



2 2 

m L m^ 



m" T m 



L"°T 



h h h 2 f h 2 

H —j{m L ) + 



m 



m T \ m L 



k A — m% 



.-, , 2 ' 1 )/,•„,./„(/,•. ////..()) • — J— p-//,, ./,,(/,-.(). //// )^ | 



2 2 

m^m L 



\N5 ab e{^[-2m 2 L j{m L ) + A; 4 (J(fc, 1, 1) + J(k,m L ,m L )) 
- (A; 2 + m 2 ) 2 J(k, m L , 0) - (A: 2 - m|) 2 J(Jfe, 0, m L ) 



A- 4 



ml 



+ 



J mn (A;, 1,1)- J m „(A;, m L , 0) - J m „(A;, 0, m L ) + J mn (A;, m L , m L ) 

k 2 ( k 2 \ k 2 ( k 2 \ 

T + 1 k m J n (k,m L ,0) j — -1 \ k m J n (k,0,m L ) 



m L \ m L 



mi \rrr L 



k 1 1 

H — T (k m J n (k,m L ,m L ) + k n I m (k,l,l)) \ >, (A. 15) 

m L J J 

n^J(A;) = ~g 2 N5 ab [u mn (k, m D , m D ) + 4k m J n (k, m D , m D ) + k m k n J(k, m D , m D )pL.16) 
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ab 



g 2 NS ab 6 mn j(m D ) 



(A.17) 



Diagrams a), b) and c) were calculated in fl9j using the Landau gauge (£ = 0). The results 
of these calculations coincide with ours if in the above formulas one sets £ = 0. The diagrams 
contributing to the 2-point function of Aq are given by diagrams f) and g): 



_ _ 




The corresponding analytical contribution is: 

k 2 + m ? 



2 N5 nb 5^5 0u 



-g ivo a6 
1 



j{m T ) 



m 



T 



m 2 D + k 2 ) 



rrin 



-J(k,m D ,0) (A.18) 



+ 



2k 2 (my + m 2 D ) + k A + (m? — m 2 D ) 2 ) J{k, tub, tut) — Hj^d) 



m. 



g 2 N5 ab 5^5^\ {m 2 L + m 2 D + k 2 ) j{m L ) + (m 2 D + k 2 ) 2 J(k, m D , 0) 



(m 2 D + k 2 ) 2 J(k,m D ,m L ) 



11 afe 



NS„ h 6 0,1 6 0v 



9 -lyOab 



2j(m T ) +ij(m L ) 



(A.19) 



Finally the single diagram contributing to the ghost 2-point function is given by diagram h): 

h 



The corresponding analytical contribution is: 

~m\ - (k 2 + m 2 G 



^ab(k) 



-\g 2 N5 ab 



+ 



(k 2 + m 2 r ) 2 T/ , 
j(m T )- y Am2 G> J(k,m G ,0) 

{k 2 + m 2 G ) 2 k 2 — m 2 G m%\ I 



Am 2 T 



(A.20) 



2""^ + — I J(k,m G ,m T ) - -j(m G ) 



(k 2 + m 2 r ) 2 k 2 + m 2 r ml \ T ,. . 1.. , 



{k 2 + m 2 G ) 2 3m 2 L + k 2 + m G 



4m| 



-J(fc,m G ,0) 



4m| 
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B 3d Adjoint Higgs Model on Lattice 



B.l Phase Diagram of 3d Adjoint Higgs Model and Dimensional 
Reduction 

As it was discussed in section 2 the high temperature limit of 4d SU(N) theory is described 
by 3d adjoint Higgs model with the following action 

S = j d 3 x^F«F« + \(D t A a ) 2 + \m 2 m A*Al + X A (A a A a ) 2 ) , (B.21) 

The relations between the parameters appearing in this equation and those of the original 



4d SU(N) theory are [|4 



where 



gl 



9 (T _L 

■2 



DO 



\{N+ l -N s )g\^)T 



1 



g- 



(6 + N — Nj 



24tt 2 



1 + 2 



16tt 2 

g 2 



16n'- 



(L + c„ 
:(L + c 



(JV)^ 



.(2) 



J3) 



22iV, a 4N fl 4u 
In in — . 

3 /Mr 3 [It 

N 
~3' 

10iV 2 + 2iV 2 + 9N f /N 

QN + 3N f ' 

7/3 - 109.N//96 2 

1 - N f /8 + 3 

7/2-23iV f /18 2 Ar 

— J - h -N f 

l-N f /9 3 f 



N f 



(B.22) 
(B.23) 
(B.24) 

(B.25) 
(B.26) 
(B.27) 
(B.28) 
(B.29) 



with \xt = 47re~ 7£; T ~ 7.0555T. One can easily verify that the parameters g%, rriDo an d 
do not depend on the renormalization scale fi at leading order. Following Ref. ]34| we 
introduce the following dimensionless parameters which will be useful in further analysis 



gi 9t 



2 

DO 



(B.30) 



The lattice action corresponding to the continium action ( B.21 ) using standard discretization 
procedure could be written as 

Sa = ^^[/p + ^2[TraAg(x)-TraA (x)C/ i (x)A (x + i)C/|(x) 



+ 



^[(m a) 2 Tra^( x ) + a\ A (Tr aA 2 ) 2 



(B.31) 
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where the first term is the standard Wilson action for the gauge fields, A Q = J2b rbj ^o with 
r b being the gararators of the SU(N) gauge group and a is the lattice spacing. Futhermore, 



P 



2N 



(B.32) 



and m is the cutoff dependent bare mass of the A field calculated in lattice regularization. 
In the following we will consider only the SU(2) case. In order to relate the results of 
lattice calculations to the physics of the high temperature phase of the original 4d SU(2) 
gauge theory, it is necessary to find the relation between the bare mass m and and the 
renormalized mass at scale \i in MS scheme m r D (fi). Since the theory is superrenormalizable 
the only running coupling is the mass parameter. The scale dependence of the renormalized 
mass is given by 

, 5(2g| -Xa)Xa " 2 - 
m D\N — m ™ H — — m • 



m D0 + 

The bare mass fho could be written as 



16tt 2 



m = m r D {y) + Sfh(fjL), 



where 



3.1759114 
4na 



(4g 2 3 + 5X A )- 



16vr 2 



(20x - 10x 2 )(ln — + 0.09) + 8.7 + 11.6s 
an 



(B.33) 



(B.34) 



(B.35) 



is the 2-loop lattice counterterm calculated in Ref. | 102 |. It is convenient to write the lattice 
action in terms of antihermitian matricies Aq defined by (3Aq = 2ia 1 ^ 2 Ao. Futhermore 
introducing the bare mass parameter h 



m^a 2 and taking into account that 



A A Tr(a^o) 2 = x/3(-Tr Aq 



2\2 



^B.36) 



the lattice action could be written in the form 

S = (3 



^2 



^Tt Up + PY, 1 Tt io(x)[/i(x)i (x + i)£//(x) + 



E 



i 



P[3 + -h) -TtA 2 (^)+(3x[ -TrAfe 



[B.37) 



Making use of Eqs. (p.34j ) and ( |B.35|) the bare mass parameter h could be written in terms 
of P and continuum parameters x and y as 



h(x) 



16 

P 



2 y 



3.1759114(4 + 5x) _ 1 _ ^ 3^ + Q _ og) + g _ ? + n x 

up ir z p z V 2 / 



The phase diagram of the 3d SU(2) adjoint Higgs model was established in Ref. p3| in terms 
of x and y using multicanonical update technique ||103|| . It has a symmetric (confinement) 
phase and the broken (Higgs) phase separated for small x by a line of 1 st order phase 
transition. The transition becomes weaker as x increases and eventually turns to a crossover 



OS 



for some x > 0.3. However, dimensional reduction is applicable only for temperatures larger 
than the critical which means that physically relevant values of x are smaller than 0.2. The 
transition line was obtained by fitting different data points with a polynomial and was found 
to be (see Figure 

xy c = 2/(9tt 2 )(1 + 9/8x + Ax 3/2 + Bx 2 + Cx 5/2 + Da; 3 ), (B.39) 

where A = 45.1(2.9), B = -214.7(20.2), C = 350.7(44.5), D = -206.7(31.4). 

The physics of the high temperature SU(2) gauge theory corresponds to some line in 
xy-plane, the line of 4d physics y^dix). This line can be calculated using Eqs.( |B.22|) -( p.24|) 
and for SU(2) it has the following form 

( 8 -N f ){A + N f ) 192-2^-7^-2^ 
MX) ~ 144A + 96(8-^ • (B - 40) 

The position of y4d(x) relative to y c (x) determines the physical phase of the 3d adjoint Higgs 
model. It turns out the for physically interesting values of x y&d{x) < Vc(x) and the physical 
phase is the broken one. The corresponding situation is shown in Figure PU[ 

The fact that the broken phase corresponds to the high temperature phase is self- 
contradictory because dimensional reduction is valid if A <C y/T, but in the broken phase 
A ~ 1/g. However, the line of 4d physics lies in the region where the symmetric phase is 



still metastable for finite lattice volume as it is indicated in Figure 20. Using this fact it was 



suggested in Ref. |33| that one can perform simulations along the line of 4d physics y^{x) 
analyze samples in the symmetric branch of the mixed phase. The obvious problem with 
this approach is that the metastable phase disappears in the infinite volume limit. 



A different approach for fixing parameters appearing in Eq. (|B.37|) was proposed in Ref. 



26fl . This approach is based on matching different quantities calculated with weak coupling 
expansion in the 4d lattice gauge theory and the corresponding 3d effective theory at some 
small physical distance, where the weak coupling expansion is applicable. In this approach 



the bare mass parameter is given by 
where Lq and L s are the temporal and spatial lattice sizes of the original 4d lattice and g is 



the 4d gauge coupling. n o was calculated in Ref. ||104|| . The other two couplings appearing 
in the lattice action are given by 

= r 2Mx 2^ (B.42) 
L g 2 (l + ag 2 ) 

(37T 2 (1 + ag 2 )) 

where a is a slightly volume dependent factor. The gauge coupling g was found by match- 
ing the Polyakov loop correlator measured in Monte-Carlo simulations to its perturbative 
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Figure 20: Data points in the limit V — > oo, a — ► for the critical curve y = y c {x) (multiplied 
by x). The thick line is the fit to the V = oo extrapolated data. The dashed line marks the 
region where the transition turns into a cross-over. The straight lines are the 4d— >3d curves 
of eq. ( |B.40|) marked by the value of Nf. The top scale shows the values of T/A^jg- (Figure 
is taken from Ref. 131 



expression at distance R = (4T) _1 [[L04 |. The numerical values of the coupling for lat- 
tices considered in Ref. p6| together with the physical temperature are summarized in the 
following Table Q 





T/T c 


9 2 


P 


h 


X 


4 x 16 3 


3.5 


1.26 


13.70 


-0.28 


0.046 


(a = -0.06) 


6.0 


1.18 


14.60 


-0.27 


0.043 


(n 00 = -0.0384) 












4 x 24 3 


2.0 


1.43 


12.25 


-0.31 


0.053 


(a = -0.059) 


3.5 


1.28 


13.54 


-0.28 


0.047 


(n 00 = -0.0431) 


6.0 


1.19 


14.48 


-0.26 


0.043 



Table 5: Parameters of the effective theories for different temperatures and lattice volumes 



according Ref. [26 



7 The numerical values of h given in pq] are not correct. The correct values of h calculated from Eq. 



(B.41) are slightly different from those given in Ref. EC 
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Now we look for the physical phase in this approach. This question was addressed in Ref. 
HI and it was concluded that the physical space is the symmetric one. In Ref. [M the 



values of h corresponding to the phase transition were determined for 4 x 12 3 , 4 x 16 3 and 
4 x 24 3 lattices and found to be h c = -0.2623(1), - 0.2803(1) and -0.30(1). Extrapolation 
to infinite volume yields —0.2943(2). Now using Eq. ( |B.38| ) we can evaluate h c from the 
result of the polynomial fit for y c (x) given by Eq. ( P-39|) at x ~ 0.047 and (3 = 13.54, 
which yields h c = —0.2733(25). Thus the values of h corresponding to the transition are 
underestimated. We can estimate the infinite volume limit of h defined by Eq. ( |B.41| ), the 



calculation yields h = —0.2699. However, Eq. ( |B.41| ) is only valid at 1-loop level, at 2-loop 
level h is given by Eq. ( P-38|) which yields —0.2789 which definitely corresponds to the 
broken phase. Therefore the conclusion that the symmetric phase is the physical reached in 
is sensitive to the approxiomation made in the analysis, and thus is unreliable. 



B.2 Extracting the Screening Masses from the Propagators 

For the analysis of the large distance behaviour of the propagators it is convenient to intro- 
duce the so-called local masses. They are defined by the following relation 

Gi{z) cosh( mi (z)(z - N z /2)) 



Gi(z + \) cosh(mi(z)(z + 1 -N Z /2)Y 

where i — D,T and the correspondig propagators are defined by Eqs. (f79f) and (j80"D and 
N z is the extension of the lattice in z direction. If the propagators show exponential decay 
starting from some value of z the corresponding local masses reach a plateau. Therefore the 
local masses are usefull in defining the starting point of the fit interval. In order to determine 
properly the fit interval correlated fit should be used together with the x 2 criteria. If one uses 
simple uncorrelated fit the corresponding values of x 2 '/d.o.f are very small which does not 
allow to choose the fit interval properly. In Figure |21| the local magnetic masses measured on 
16 2 x 64 lattice at /3 = 8, x = 0.09 and y = 0.4007 are shown. The filled squares denotes the 
data points which were used in the fit and yields good x 2 /d.o.f '. The fit based on this data 
points yields the value 0.422(6) (^f for the magnetic mass. In Ref. |9]| the magnetic mass 



was measured in the symmetric phase of the 3d SU(2) Higgs model. The magnetic mass for 
the 3d pure gauge theory was also determined there and its value was found to be somewhat 
larger but close to the corresponding value of the 3d Higgs model. In Ref. |)3| the magnetic 
masses were extarcted from propagators using uncorrelated fit over the interval [8, N z — 8] 
with the following fit anzatz 

G T (z) = A\exp(-mrz) + exp(-m r (iV 2 - z))J + B (B.45) 

where the propagators were measured for (3 = 9 on 16 2 x N z lattice, with N z = 32 — 128. For 
N z = 64 the value of the magnetic mass was found to be 0.39(2) #f (see Table 1 in Ref. ||93|| ). 
Using the same procedure for our data we obtaine = 0.408(8)^1 which is compatible 
with the result of Ref.|93|. However, most of simulations in Ref. |53] were done on 16 2 x 32 
lattice. On this lattice the magnetic mass was found to be 0.35(l)p§. 

In order to get in contact with the results of Ref. [ 93| we have performed simulations on 
16 2 x 32 lattice at (3 = 9 in 3d pure gauge theory and extarcted the magnetic mass using 
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Figure 21: The local magnetic masses in units of g\ as function of z measured at /3 — 8, a; = 
0.09 and y = 0.4007. The filled squares denote the data points used for extracting the 
magnetic screening mass 

the procedure outlined above. The corresponding value of the magnetic mass was found to 
be 0.39(3)gf. This value is smaller than corresponding value obatined by us for 16 2 x 64 
lattice and this is probably due to the fact that on such lattices the magnetic mass cannot 
be reliably determined because the local masses do not show a plateau (see Figure 2 in Ref. 
p3|). We think that the small difference between this value and the corresponding value 



found in Ref.|93| is probably due to the presence of the Higgs fields. 



The above analysis shows the importance of the finite size effects and the proper procedure 
of extracting the screening masses. 

B.3 The Dependence of the Screening Masses on the Choice of 
the Temperature Scale and Lattice Spacing 

As it was already explained in section 4.1 the temperature scale in the effective theory is set 
by parameter x and its relation to the renormalized 4d gauge coupling g(T). In section 4.1 
g(T) was determined by the 2- loop formula in MS scheme with fi = 2ttT. It seems to be 
interesting to examine to what extent the numerical results on the screening masses presented 
in section 4.3 depend on the definition of g(T) and/or the choice of the renormalization scale. 
In Figure [22] the temperature dependence of the screening masses is shown for the case when 
the temperature scale is set by the 1-loop level renormalized coupling 

S- 2(T) = l2^ ln X= (B.46) 

with jj, = 2ttT and = T c /1.06. As one can see from Figure [22] the screening masses are 
systemetically overestimated compared to the 4d results in this case. 

It was mentioned in section 4.4 that the square root of the spatial string tension in the 
deconfined phase of the 4d SU(2) gauge theory is expected to scale with g\ ~ g(T) 2 T. This 
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Figure 22: The screening masses in units of the temperature in the case when the temperature 
scale is set by the 1-loop level renormalized coupling (see text for details) Shown are the 
Debye mass mp for the first (filled circles) and the second (open circles) set of h, and the 
magnetic mass for the first (filled squares) and the second (open squares) set of h. 
Also shown there are the values of the Debye mass measured in metastable phase along the 
perturbative line of 4d physics (open triangles). The line (a) and line (b) represent the fit 
for the temperature dependence of the Debye and the magnetic mass from 4d simulations 
from |49| . Some data points at the temperature T ~ 260T C and ~ 30000T C have been shifted 
in the temperature scale for better visualization. 



fact can be used for non-perturbative definition of g(T). In Ref. 91 it was found that the 
temperature dependence of the square root of the spatial string tension is well described by 
the formula 

\/a s (T) = (0.369 ± 0.0U)g 2 (T)T (B.47) 
if the running coupling constant g is chosen according to 



9-\T) 



11 T 17 
In 1 

12vr 2 A T 44tt 2 



In 



2 In 



T 



;B.48) 



with At = 0.076(13)T C |9T|. This running coupling constant g(T) is equivalent to the 2-loop 
running coupling in MS scheme if the renurmalization scale is chosen to be fi ~ 13. 3T. 
However the data on the temperature dependence of the spatial string tension can be also 
well fitted with the formula ^Ja s (T) = (0.334 ± 0.014)g 2 (T)T if g(T) is chosen according to 



9-\T) 



11 T 
?r In -; — 

127T 2 At 



(B.49) 



with At = 0.050(10)T C ||91|| ■ This running coupling corresponds to the running coupling 
constant in MS scheme with \x ~ 20T. The value of the string tension in 3d SU(2) gauge 
theory is ^fal = 0.3340(25)^1 [p2| . Now if the temperature scale is set by the running 
coupling g(T) defined by Eq. (|B.46|) but with /x ~ 20T very good agreement between the 3d 



and 4d data is found. 

From the above analysis it is clear that the definitions of g(T) and values of the nor- 
malization scale /i which yield good agreement between 3d and 4d data are compatible with 
values of g found non-perturbatively from the analysis of the spatial string tension. On the 
other hand the values of gauge coupling used in Ref. |EB| are by factor of 2 smaller than 
these. 
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The matching analysis in section 4.3 was done at (3 — 16. It is important to check 
the sensitivity of these results to the choice of /3 which defines the lattice spacing. The 
dependence of the Debye masses on (3 was analyzed for x = 0.05 and for all three values of 
y (h) used in our matching analysis. The corresponding results are shown in Table 6. 



p 


12 


16 


20 


y = 0.6721 


1.84(21) 


1.69(6) 


1.68(12) 


y = 0.5914 


1.57(6) 


1.68(4) 


1.65(6) 


y = 0.4756 




1.48(5) 


1.44(12) 



Table 6: Debye masses for different values of (3 for x = 0.05 

As one can see from the results presented in Table 6 no lattice spacing dependence for the 
Debye mass can be observed. A similar analysis should be done for the magnetic mass. But 
for (3 = 20 the magnetic mass can be measured reliably only for lattices larger than 32 x 64. 
Using the fact that the magnetic mass calculated in the 3d adjoint Higgs model is close to 
the magnetic mass of 3d pure gauge theory one and based on results presented in section 4.4 
we can conclude that the magnetic mass also does not depend on (3. 

C Self-energy Contributions in 3d fundamental + ad- 
joint Higgs model 

Below we list graphycally the additional diagrams contributing to the A (a-i), Higgs boson 
(j-1) and the vector boson (m) self-energies and the vacuum expectation value (n). Here solid 
lines denotes the propagator of the Higgs field, dashed lines corresponds to Goldstone fields, 
double line corresponds to the adjoint Higgs field and finally the dotted line corresponds to 
ghost fields. Wawy lines as usually corresponds to vector fields. 




a be 




d e f g h i 




m n 
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